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Abstract: We consider control-constrained linear-quadratic optimal control problems on evolving 
hypersurfaces in In order to formulate well-posed problems, we prove existence and uniqueness 

of weak solutions for the state equation, in the sense of vector-valued distributions. We then carry 
out and prove convergence of the variational discretization of a distributed optimal control problem. 
In the process, we investigate the convergence of a fully discrete approximation of the state equation, 
and obtain optimal orders of convergence under weak regularity assumptions. We conclude with a 
numerical example. 
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1 Introduction 

We investigate parabolic optimal control problems on evolving material hypersurfaces in 
R n+1 . Following [DE07], we consider a parabolic state equation in its weak form 

S j t j y^&T(t)+ j Vri/-Vr¥>dr(t)= j y^dT(t) + j fip6T(t), (1.1) 



o 



r(t) r(t) r(t) r(t) 



where T = {r(t)} L ' J is a family of C 2 -smooth, compact n-dimensional surfaces in R n+1 , 
evolving smoothly in time with velocity V. Further assume / sufficiently smooth and let 
ip = dt<p + Wip denote the material derivative of a smooth test function if. 
We start by defining unique weak solutions for the state equation. The idea is to pull back 
the problem onto a fixed domain, introducing distributional material derivatives in the sense 
of [LM68] and a W(0, T)-like solution space. As a consequence, a large part of the theory 
developed around W(0, T) for fixed domains applies, compare for example [LM68] and [Lio71] 

An alternative approach to prove existence of weak solutions along the lines of [LSU68] is 
taken in [SchlO], that entirely avoids the notion of vector- valued distributions. 
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Recent works also deal with the discretization of (1.1), both in space, compare [DE10], and 
time, see [DLM11] and [DE11]. 

In [DE10] order-optimal error bounds of type sup te r 0T i II ' llL 2 (r(t)) are derived for the dis- 
cretization of the state equation, assuming a slightly higher regularity of the state than is 

used in section 5 and 6, where we derive ^ J Q T || • H^p^dt^ 5 -like bounds. A class of Runge- 

Kutta methods to tackle the space-discretized problem is investigated in [DLM11], assuming 
among other things that one can evaluate / in a point-wise fashion, i.e. that f(t) G L 2 (T(t)) 
is well defined. For a fully discrete approach and the according error bounds see [DE11]. 
There a backwards Euler method is considered for time discretization whose implementation 
resembles our discontinuous Galerkin approach in Section 6. Yet while the approach in [DE11] 
ultimately leads to sup t6 [ T ] || • || L 2(- r ( t ))-convergence, we allow for non-smooth controls and 
thus cannot expect to obtain such strong convergence estimates. 

Basic facts on control constrained parabolic optimal control problems and their discretization 
can be found for example in [Tr605] and [MV08b], respectively. 

The paper is structured as follows. We begin with a very short introduction into the set- 
ting in Section 2. In order to formulate well posed optimal control problems we first proof 
the existence of an appropriate weak solution in Section 3, complementing the existence re- 
sults from [DE07]. We then use the the results from Section 3 in order to formulate control 
constrained optimal control problems in section 4. Afterwards, we examine the space- and 
time-discretization of the state equation in Sections 5 and 6, before returning to the opti- 
mal control problems in Section 7. There we apply variational discretization in the sense of 
[Hin05] to achieve fully implementable optimization algorithms. We end the paper by giving 
a numerical example in Section 8. 



2 Setting 

Before we can properly formulate (1.1), let us introduce some basic tools and clarify what 
our assumptions are regarding the family {r(*)} t6 [ T]' 

Assumption 2.1. The hypersurface To = T(0) C R n+1 is C 2 -smooth and compact (i.e. 
without boundary). F evolves along a C 2 -smooth velocity field V : M n+1 x [0, T] — > W l+1 
with flow $ : R n+1 x [0,T] 2 ->• R™ +1 , such that its restriction $£(•) : T(s) ->• T(t) is a 
diffeomorphism for every s,t G [0, T\. 

The assumption gives rise to a second representation of F(t) and in particular implies T(t) 
to be orientable with a smooth unit normal field u(-, t). As a consequence, the evolution of V 
can be described as the level set of the signed distance function d such that 

r(t) = {x G R n+1 I d(x,t) = 0} , 

as well as |d(ar, t)| = dist(x, T(t)) and Vd(x,t) = u(x,t) for x G r(£). Further, we have 

d(-,t) G C 2 {M r {t)) for some tubular neighborhood M r (t) = {x G R n+1 \ \d(x,t)\ < r) of Y{t). 

Due to the uniform boundedness of the curvature of T(t) the radius r > does not depend 

on t G [0,T]. The domain of d is Mt = Utg[o T] N r (t) x {t} which is a neighborhood of 

U e[ o,T]r(t)x{t}inR™+ 2 . 

Using d we can define the projection 

at : Kit) ->• T(t) , a t (x) = x - d(x,t)Vd(x,t) , (2.1) 
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which allows us to extend any function <f> : F(t) — > R to Af r (t) by 4>{x) = 4>{at(x)). Hence we 
can represent the surface gradient in global exterior coordinates Vm)^ = {I ,t)v{- } t) T )\/ cf) 
as the euclidean projection of the gradient of 4> onto the tangential space of T{t). In the fol- 
lowing we will write Vr instead of Vr(t), wherever it is clear which surface Y{t) the gradient 
relates to. 

We are going to exploit existing results on vector- valued distributions, which we recall 
here for completeness. In order to define weak derivatives consider V((0,T)), the space 
of real valued C°°-smooth functions with compact support in (0, T). Fix s € [0, T]. Each 
y e L 2 ((0,T),H 1 (T(s))) defines a vector-valued distribution T y : V((0,T)) -»• i? 1 (r(s)) 
through the ff 1 (r(s))-valued integral J y(t)ip(t)dt. 

[o,n 

Its distributional derivative is said to lie in L 2 ((0, T), H^ 1 (r(s))) iff it can be represented by 
w e L 2 ((0,T),H^ 1 (T(s))) in the following sense 

7»= j y(t)<p'(t)dt = - J witMQdteH 1 ^)), V<p€V((0,T)), (2.2) 

[0,T] [0,T] 

and we write y' = w. Note that by R~ x we denote the representation of the dual (H 1 )* which 
arises from L 2 D H 1 by completion. 

Lemma 2.2. For s € [0, T], i/te space 

w s (o,r) = {d£ ^((crj.^HrW)) | «' e ^((o.r),^- 1 ^)))} 

loit/i scalar product Jq (■, •)// 1 (r(s)) + ((")'> (')')j : /- 1 (r(s))dt a Hilbert space. 

1. W s (0,T) is compactly embedded into C([0,T], L 2 (T(s))), the space of continuous L 2 - 
valued functions. 

2. Denote by V([0,T}, H l (T(s))) the space of C°° -smooth H l {T(s)) -valued test functions 
on [0,T]. The inclusion V{[0,T], i? 1 (r(s))) C W s (0,T) is dense. 

3. For two functions v,w € W s (0,T) the product {v(t),w(t)) L 2^ r ^ is absolutely continu- 
ous with respect to t £ [0, T] and 

d f 

— / v(t)w(t) dr(s) = (v, w) H -i ir{s))tH i {r{s)) + (u, ^%i(r( s )),H-i(r( s )) , 

a.e. in (0,T), and as a consequence there holds the formula of partial integration 
J (v',w) H -i tH i dr = (v(t),w(t)) L 2 (r{s)) - (v(r),w(r)) L 2 {r(s)) - J {v, w') m>H -i dr . 

[r,t] [r,t] 

For a proof of the lemma, see [LM68, Ch. I,Thrms. 3.1 and 2.1]. In fact one can use the 
formula of partial integration to prove the embedding into C((0, T), L 2 (T(s))), see [Eva98, 
Ch. 5,Thrm 3]. For further references see [TV605, Thm. 3.10]. 



3 



Our approach to weak material derivatives relies on the following equivalent formulation of 
condition (2.2) 



J (y(*),¥' / (*))L2 ( r( s )) + (^(t),¥'(t))ff-i(r( s )),ffi(r( s ))dt = 0, ( 2 - 3 ) 

[0,T] 

which defines the weak derivative y' = w of a function y G L 2 ((0, T), // 1 (r(s))) via its 
L 2 ((0,T),L 2 (r(s)))-scalar product with elements of D((0, T), i7 1 (r(s))). 
The equality (2.3) follows from (2.2) by Lemma 2.2 [5.]. On the other hand (2.2) is a conse- 
quence of (2.3). To see this, test (2.3) with G Z>((0, T), i? 1 (r(s))), where V G X>((0,T)) 
and v G fl" 1 ^))- 



3 Weak solutions 



The scope of this section is to formulate appropriate function spaces and a related weak 
material derivative, in order to prove the existence of unique weak solutions of (1.1) for quite 
weak right-hand sides /. 

We start by defining the strong material derivative for smooth functions / G C 1 (R n+1 x [0, T]), 
namely the derivative 

/(z,t) = -f f(<S>l(x),s) = Vf(x,t)V(x,t) + d t f(x,t), (3.1) 

US s=t 

along trajectories of the velocity field V. The material derivative has the following properties. 
Lemma 3.1. Let f be sufficiently smooth. Then 

d 
dt 



J /dT(t) = J /+ divr^dr(t), 



r(t) r(t) 

and 

d_ 
dt 



J ||V r /|| 2 dr(t) = J 2Vr/-Vr/-2Vr/( J Dr^)Vr/+||V r /|| 2 div r ydr(t) 



r(t) r(t) 
with div r{t) V = Er=i' V r(t) ^ and (D r(t) V) t] = V J r{t) V\ 
A proof and details can be found in the Appendix of [DE07]. 

Lemma and Definition 3.2. Let Jf{ ■ ) = detD r / s )$|( • ) denote the Jacobian determinant 
of the matrix representation o/ Dp( s ) <&f ( ■ ) with respect to orthogonal bases of the respective 
tangent space. By Assumption 2.1 Jf G C 1 ([0, T] x T(s)) and there exists Cj > 0, such that 
for all s, t G [0, T] 

±- < min J/( 7 ) < max J/( 7 ) < Cj . 
Cj 7Gr(s) 7er(s) 

Given Assumption 2.1, consider the family {L 2 (r(t))} te j Q ^ . Then for v € L 2 (T(t)) we in- 
troduce the pull-back 

4> s t v = v(<f>i(-))eL 2 (r(s)), 



4 



which is a linear homeomorphism from L 2 (T(t)) into L 2 (F(s)) for any s,t G [0, T]. Moreover 
(f>f is a linear homeomorphism from // 1 (r(t)) into // 1 (r(s)). 77ms finally the adjoint oper- 
ator, (j)** : H~ 1 (T(s)) — > i7 _1 (r(i)) is also a linear homeomorphism. There exist constants 
C L 2(i>j, C^i(r) independent of s,t, such that for allv G L 2 (T(t)), or v G // 1 (r(t)) respectively, 
and for all s, t G [0, T] 

II^H/PfTOO) ^ CtfHnlMI-ffHrX*)) ' ll$ u IU 2 (r(s)) < <^L 2 (r)IMlL 2 (r(t)) , 

and /maZZt/ \\(f>i*\\c(H-i(r(s)),H-i(T(t))) < Cffi(r)- 
Furthermore there holds dtJf = (/>*( div r(t)^)<^t • 

Proof. For s, i G [0, T] we have 

| « 2 dr(t)= | (^) 2 J/dr( s ) 

r(t) r( s ) 

i 

and thus 1 1 0f ^ II z.2 ( r( s )) < C L *(r) ||«||L»(r(t))' with ^(r) = Cj- 

For ^f 1 equivalence consider v G i? 1 (r(t)) and choose G C 1 (r(s)). Now 

1 (#t;)Vr¥>dr( s )= 1 t;(D^) T Vr(^^)J*dr(t) 
r(s) r(t) 

and because -u G i^ 1 (r(t)) we can integrate by parts on T(t) to obtain with u s = s) 

y (# V )Vr¥>dT( a ) = - y «; (^^)J S * dT(t) = - y - H s v s &v)v + {4> s t v)^H s u s dT(s) . 
v{s) r(t) r( s ) 

Note that u; G L 2 (r(t)) n+1 and that ||«;||x,2rr(t))«+ 1 — CIMIffifiYt))) where C depends only 
on the mean curvature H t of T(t) and the second space derivatives of $ which are bounded 
independently of s,t G [0,T]. Now Vr( s )(</>£^) = — H s v s 4>\v G L 2 (r(s)), because as stated 
above ||</>HlL 2 (r(s))™+i < C L 2 {r) C\\v\\ H i {T{t)) , and || L 2 (r(s)) ^ C L 2 (r) |M| L _ 2(rW) . Thus, for 
some C H i(r) > depending only on a global bound on \H t \, ||5j$f|| and ||9jj$|||, s,t G [0, T], 
1 < < 72 + 1, there holds 

1 1 ^11 //MrM) ^ c m(r)\\v\\m(v(t)) ■ 

Now || • ||/fi(r(t)) an d 110? ( " )llH 1 (r(s)) are two equivalent norms on i? 1 (r(t)). Hence also their 
dual norms are equivalent. The norm of / G (i^ 1 (r(s)))' can now be expressed by 

Qn „ (f> w )(m(r(s))y,m(r(s)) _ (<l>t*f,v) HHm) /« 9 x 

sup — sup , ^.ZJ 

men 1 ^*)) lFllHi(r( s )) «eHi(r(t)) ll0t«l|Hi(r(s)) 

and the bound on the norm of (f)f * follows from the equivalence of said TZ^-norms. 

The last assertion is a by-product of the proof of Lemma 3.1, compare [DE07]. I 

We need to state one more Lemma concerning continuous time-dependence of the previously 
defined norms. 
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Lemma 3.3. Let s G [0,T]. For vi G i? 1 (r(s)), v 2 G L 2 (T(s)), v 3 G fr _1 (r(s)) the following 
expressions are continuous with respect to t G [0, T] 

H^illtfHrxt)) , ||$W2||i/«(r(t))> ll#* u 3||jr-i(r(t))- 
Proof. By the change of variables formula we have 

ll#«iHW( t ))= / (VmPrw^TH^^r'Vm + ^^drM, (3.3) 
r( s ) 

which is a continuous function due to the regularity of <1> stated in Assumption 2.1. Similarly 
we conclude the continuity of the L 2 -norm. 

Moreover, since (D^ir^D^^ = id 7T( s ), ^ = 1, and G C 2 (r(s) x [0,T]) 

Equation (3.3) infers 

lll^s^lllf^r^)) ~ IMIff!(r(s))l - ^1* ~ s lll' u llff 1 (r(s)) ' 
for all v G i^ 1 (r(s)). Regarding (3.2) this allows us to estimate 

-—^ — -7TlNlff-i(r( a )) ^ ll^r«3||H-i(r(t)) < - — — — t\\^\\h-\y{s)) ■ 
(l + C\s — t\)2 (1 — C7|s — t|) 2 



As far as Lemma 3.1 is concerned, for a family of functions {/(i)}te[0,T]> /(*) • — >■ M, one 
can define / at 7 = $070 simply by /(t)[ 7 ] = $&($/(*))l7o,*] = $£[/(*)(* t °7o)]. If {/(*)} 
can be smoothly extended, this is equivalent to (3.1). The following Lemmas aim at defining 
a weak material derivative of / that translates into a weak derivative of the pull-back (f>tf(t). 

Lemma and Definition 3.4. Consider the disjoint union B L i = [j t ^ 0T ^ L 2 (F(t)) x {t}. 
The set of functions f : [0,T] — > B L 2, t i-> (v t ,t) inherits a canonical vector space structure 
from the spaces L 2 (T(t)) (addition and multiplications with scalars). Given Assumption 2.1, 
for s G [0, T] we define 

L h(F) ■= {v:[0,T]^B L 2,t^(v t ,t) \ <f>\v G L 2 ((0, T), L 2 (T(s))) } . 

Abusing notation, now and in the following we identify v(t) = (vt,t) G L 2 ^^ with v(t) = Vt- 
Endowed with the scalar product 

</,<?>L2 2(r) = J (f(t),g(t)) LHm) dt. 

[0,T] 

L 2 2 ^ r -j becomes a Hilbert space. 

In the same manner we define the space L 2 Hl ^ r y For L 2 H _ 1 ^ use (f) 1 * instead of <fif. All three 
spaces do not depend on s. 

For ip G 4' ) P((0,T),F 1 (r(s))) = {<p G L 2 2(r) | <ft<p G D((0, T), H 1 ^)) }, it is clear how 
to interpret (p, namely (p = fiKfifip)' G i7 1 (r(t)). We say that y G L 2 Hl ^ has weak material 
derivative y(t) G L 2 H _ 1 ^ iff there holds 

J <J/,¥>)ff-i(r(t)),tf 1 (r(t)) d * = ~ / (y'^>L 2 (r(t))dt- y ^ div r V^dr(t) dt (3.4) 

[0,71 [0,71 [0,T] r(t) 
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for all (p G (j)i^T>((0,T),H 1 (T(s))), and the definition does not depend on s. 

Proof. In order to define the scalar product of L 2 L2 ^, we must ensure measurability of 

</,<7>L'(r(f)) : [°> T ] -> M - Since (/>0> = UWf + 9\\ 2 ~ ll/ll 2 " llsf) it suffices to show 
measurability of ||/||£2(r(t)) f° r an / e ^i 2 (rv ^ definition of the set L 2 L2 ^ we have 
0|/ G L 2 ([0, T], L 2 (r(s))). Hence, there exists a sequence of measurable simple functions 
f n that converge pointwise a.e. to <\>\f in L 2 (T(s)). Each f n is the finite sum of measur- 
able single-valued functions, i.e. /„ = fi, n l Bi , M n G N, / i>n G L 2 (r(s)), [0,7] D B< 
measurable and disjoint. By Lemma 3.3 the function 

M n 

Il0l/«||z,2(r(t)) = H^s/i,nlU 2 (r(t)) 1 B i 
i=l 

is the finite sum of measurable functions and thus measurable. Using the continuity of the 
operator as stated in Lemma 3.2, one infers pointwise convergence a.e. of ||0|/«||x, 2 (r(t)) 
towards ||/||L 2 (r(t)) which in turn implies measurability of ||/||i,2(r( t )). 

Again by Lemma 3.2 we now conclude integrability of ||/||i,2( r ( t ^ and at the same time 
equivalence of the norms 



\ 1 / \ - 

\ 2 / \ 2 



2 At 

L 2 (r(*)) m 



and j / ||#/||£ 2(r(i)) d* 
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Completeness of L 2 L2 ^ follows, since L 2 L2 ^ and L 2 ((0,T), L 2 (T(s))) are isomorph. Again 
because of Lemma 3.2, <\>%v G L 2 ((0, T), L 2 (r(s))) is equivalent to <f>\v G L 2 ((0, T), L 2 (T(r))), 
thus the definition does not depend on the choice of s. For L 2 R1 ^ and L 2 H _ 1 ^ we proceed 
similarly. 

We show that the definition of the weak material derivative does not depend on s G [0,7]. 
On r(s) Equation (3.4) reads 

J {<fc*V, <P)H-i(r(s)),m(r(s)) dt = - J J (<fty<ff '(t) + # (y div r(t) V) (p) J? dT(a) dt 

[0,31 [0,T] T(s) 

(3.5) 

for all (p G T^([0, T], 77 1 (r(s))). For r G [0, T], we now transform the relation into one on T(r), 
using ft, (<j> s r )* and <% = <j) r s o $ 

| <<#*y, ^0)H-i(r(r)),ffi(r(r)) dt = " / / ^ + ^ ( y div r « F ) ^ dI » dt ' 

[0,T] [0,T] r(r) 

and because <j) T s : 77 1 (r(s)) — > i7 1 (r(r)) is a linear homeomorphism, it also defines an isomor- 
phism between V([0,T\,H 1 (T(s))) and V([0, T], i7 1 (r(r))). I 

Remark 3.5. Strictly speaking the elements of L 2 X ^ are equivalence classes of functions 
coinciding a.e. in [0,T], just like the elements of L 2 ((0, T),X(F(s))). 

The definition of the weak derivative of y G L 2 R1 ^ in (3.4) translates into weak derivatives 
of the pullback (f>fy. In order to make the connection between the two, we state the following 
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Lemma 3.6. Let w G W s (0, T) and f G C 1 ([0,T] x T(a)). Then fw also lies in W s (0,T) and 

(fw)' = dtfw +fw' , 

6L2([0,T],L2(r( S ))) 

ui/iere /u/ is to 6e understood as (fw 1 , <p) H -i(r(s)),m(r( s )) = ( w ' J<f) H-i(r(s)),m(r( s ))- 

Proof. We show that for ip G P((0, T), i7 1 (r(s))) the function lies in W s (0, T). The claim 
then follows by partial integration in W s (0,T). 

1. Because / G C([0, T] xT(a)) and the strong surface gradient V r(s) / G (C([0, T] x r(s))) n+1 
are continuous and thus uniformly continuous on the compact set [0, T] x we infer 
/ G C([0,T],C 1 (r(s))). Note that dist [0iT]xr(s) ((t, 7), (t + fc, 7)) = k. Let e > 0, then for 
sufficiently small k e > k > one has 

n+1 

||/(t + k, •) - f(t, -)||oo + £ H V ^)/(* + M " V J r(s) /(t, OHoo < e . 

i=l 

2. As to the distributional derivative of /<£>, we show that / G C 1 ([0, T], C(T(s))). Observe 
that the uniform continuity of the strong derivative dtf on [0, T] x F(s) allows us to estimate 

\\f(t + k, •) - /(*, •) - $/(t, OfcHoo = || A: J d t f(t + rk, •) - $/(t, •) drlU < ke 

[0,1] 

for k e > k > sufficiently small. Again by uniform continuity of dtf we conclude dtf G 
C([0, T], C(r(s))). All told, taking into account the continuity of the pointwise multiplication 
between the respective spaces, we showed 

f v g c([o, T\,HHr(s))) n ^([o, t], L 2 (r( s ))) c w 9 (o, t) . 

3. Consider now an arbitrary w G W s (0,T). Since ftp G W s (0,T), by partial integration as 
in Lemma 2.2 [3.] it follows 



J («''.MH-i(r( > )),fli(r( > ))dt = - y (w, (/^)')ifi(r( s )),ff-i(r( s )) d * 

0,T] 

y (w,d t f<p) L 2(r(s))dt- y («>, fv')i?(T{s)) d * ■ 



[0,T] [0,T] 



[0,T] [0,T] 

Reordering gives 

y (fw,(p') L 2 {r(s)) dt = - J (dtfw + fw',ip) H -i (r(s))tH i (r{s)) dt 

[0,T] [0,21 

for any y> G P((0, T), // 1 (r(s))) . Hence condition (2.3) holds for fw. Using the density 
property stated in Lemma 2.2[2.], we can approximate fw by continuous i/ 1 (r(s))-valued 
functions and infer fw G L 2 ((0, T), i7 1 (r(s))). The same argument yields dtfw + fw' G 
L*((0,T),H-i(T( S ))). ' ^ I 
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Finally we can define our solution space. 

Lemma and Definition 3.7. The solution space Wr is defined as follows 

W r = {v€ L 2 Hl{r) | v G L|_ 1(r) } . 

Wr is Hilbert with the canonical scalar product Jq{-, •)_f/ 1 (r(t)) + ((")> (')) //- 1 (r(t))di- Also y G 
Wr iffffiy G W s (0,T) for (every) s G [0,T]. For a// G T>((0, T),H 1 (T(s))) there holds 

J (<f> t s *y,<P)H-i(r(s)),Hi(r(s))dt= J {((<f>ty)' ,J?<p)H-i(T(s)),Hi(T(s))dt- (3-6) 

[0,T] [0,T] 

One /ias 

cwII0?2/IIvk s (o,t) < ||y||w r < CwMy\\Ws(0,T) , 
and cw,Cw > do not depend on s G [0,T]. 

Proof. For y G W r , observe that Jf<pfy G F 2 ([0, T], i7 1 (r(s)) and rewrite (3.5) as 
J {J'tiVidtflvMs)) dt = ~ J (&*V'V)H-\r{8)),H\T{ay)te- J ( d tJt s ftty,&mr(s)) d *, 

[0,T] [0,T] [0,21 

(3.7) 

for G f ((0, T), i? 1 (r(s))). Hence Jt<frtV ^ W s (0, T), and from Lemma 3.6 it follows that also 
4> s t y G W s (0,T), because ^ G C l ([0,T] x T(s)). Note that we used d t J t s = 4> s t { dw T (t)V)J t s , 
see Lemma 3.2. On the other hand, for any y G W s (0,T) one has Jfy G W s (0, T) and thus 
y = fisV ^ Wr- Hence </>^ constitutes an isomorphism between Wr and W s (0,T). 
Apply Lemma 3.6 a second time to obtain ( J^)' = dfJfip + and because of </3(0) = 
</3(T) = G i^ 1 (r(s)) by partial integration there follows from (3.7) 

J (<l> t s *y,<p)H-i{r(s)),m(r(s))<tt= J <((#J/)' > J tV)i2-i(r( s )),ifi(r( s )) dt , 

[0,T] [0,T] 

compare Lemma 2.2 [5.]. This proves the second claim. 

The claim of Wr being Hilbert now follows. Observe that point-wise multiplication with Jf 
constitutes a linear homeomorphism in H l (F(s)) whose inverse is the multiplication by -jg. 

One easily checks \\Jt<p\\w(r(s)) < c ll^# lloi(r(s)) II II if i (r( s )) < C|Ml.ffi(r( s ))- This together 
with Lemma 3.2 yields the equivalence of the two norms on Wr 



/ 



MlWw) + \\y\\ 2 H-Hr(t)) dt and J \\<Pty\\Hi(r(s)) + \Myy\\H-i(r(s)) dt - 

[0,T] [0,T] 

Completeness of W s (0, T) then implies completeness of Wr- I 



Remark 3.8. Formula (3.6) can be seen as a generalization of the following relation. Assume 
<t> s t yeV{{Q,T),H\T{s))). Then 

J (4> t s*y,¥)H-^(T(s)),m(T{ s ))&t= j (y, <^)i 2 (r(t)) = ^ J t V) L 2 (r(s)) dt . 

[0,21 [0,T] [0,71 
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Using Lemma 3.6 and 2.2, it is now easy to proof 

Lemma 3.9. For two functions v,w G Wr the expression (v(t),w(t)) L 2^ r ^ is absolutely 
continuous with respect to t G [0, T] and 



J vwdT(t) =(v,w) H -i {T(t))jH i im) + ... 

(v,w) H i {r(t))jH -i {r{t)) + J vw div r(t )Vdr(t), 



dt 

r(t) 



r(t) 

a.e. in (0, T), and f/iere ao/ds the formula of partial integration 



J {v, w) h -i(t(t)),h1(t(t)) dr =(v, w) L 2 {r(t)) - (v, w) L 2 (r(s)) . . . 

J (v,«') JT i(r( T )) ) jf-i(r(T))+ y div r Vdr(r) 



[s,t] 

1 dr. 

[«,t] T(r) 



We can now formulate (1.1) in a weak and slightly generalized manner. Let b G C 1 ([0, T] x To) 
and b = 4>yb. We look for solutions u G Wr that satisfy y(0) = yo G ^ 2 (Lq) and for / G L^_ 1( - r ^ 



^ y yc^dr(t)+ y V r y V r ^ + 6y^dr(t) = (^,y) if -i (r(t)))Jf i (r(t)) 

r(t) r(t) ( 3 - 8 ) 

+(/> <P) H-t(r(t)),m{T(t)) > 



for all 99 G Wr and a.e. t G (0, T). One may equivalently write (3.8) as 

y + A m y + y ( div r(t) F + 6) = / in H-\T(t)) 

for a.e. t G (0, T). We apply known existence and uniqueness results for the pulled-back 
equation to prove 

Theorem 3.10. Let f G £#-i/ r ); 2/0 G ^ 2 (ro)- There exists a unique y G Wr, suc/t i/iai 
(3.8) is fulfilled for all (ft G Wr and a.e. t G (0, T). There holds 

\\y\\w T < C (jMlL2(r ) + Il/H^_ 1(r) ) • 
Proof. Let us relate equation (3.8) to the fixed domain T(s) via 

^ y ^J/dI\ S ) + y (Vrypr( s )^)" 1 (^r( s )^)" r Vr^ + 6^) J t s dr( s )... 
r(-) r( s ) 

= Jty)H-i(r(s)),m(r(s)) + (f, Jt <P)H-i(r(s)),m(r(s)) , 
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with y = (f)fy, f = jscpl*/ G L 2 ((0, T), H~ 1 (T(s)) and for all <ft<p = (p G W s (0, T). This again 
is equivalent to 

(y, <pJ?)H-i(r(s)),Hi(r(s)) + J VP (^(div r( t)F) + &) J/ dl» + . . . 

T(s) 

+ y V r y( J Dr(s)^)" 1 ( J Dr (s )^)" T Vr^^dr( s ) = </, J/$ff-i(r W)rff i(r W ) • 
r( s ) 

With V = one gets for all V G W a (0, T) 

W, V') J ff-i(r(s)), J ffi(r(s)) + y, V 7 ) = (/, ^>H-i(r( s )),Hi(r( s )) , ( 3 - 9 ) 
with a bilinear form 

a(t,y,i;)= J V T y(Dr^i)~\Dr{s)^i)~ T ^ dT(s) + | j/ (^(div r(t) K) + &) ^dT(s) ... 
r( s ) r(s) 

r( s ) 

By Assumption 2.1 the bilinear form (D r ^ s ^f)~ 1 ['y}(D-p( s )^t)~ [7] is positive definite on the 
tangential space T 7 r(s) uniformly in s,t G [0, T] and 7 G T(s). Thus, there exists c > 
such that for some fco > one has a(t,tp,i/j) + 1 1 V 7 1 1 z, 2 (r(s)) > c IIV ; l|j : f" 1 (r(s))- We are now in 
the situation to apply for example [Lio71, Ch. Ill, Thrm. 1.2], to obtain a unique solution 
y G W s (0,T) to equation (3.9) for initial data (fr^yo G L 2 (T(s)). Moreover the solution map 
is continuous 

||y||w s (o,T) < C (||/||L2(( 0jT ) jff -i( r ( s ))) + ||0oyo||L 2 (r( s ))) 

Note again that ||/||L2((o,T),_ff- 1 (r(s))) — C||/IIl 2 1 > smce the multiplication with Jf is a 

globally bounded linear homeomorphism in i7 1 (r(s)), as stated in the proof os Lemma 3.7. 
The transformation of (3.8) into (3.9) works both ways, hence the uniqueness of y G Wr- The 
norms can be estimated as in Lemma 3.2 and Lemma 3.7 and the theorem follows. I 

With regard to order-optimal convergence estimates, sometimes a slightly higher regularity 
than y G Wr is required. Assuming / G L 2 L2 ^ and j/o G ^(Tq), one can apply a Galerkin 
approximation argument, see [DE07, Thms. 4.4 and 4.5] for manifolds or [Eva98] for open 
sets, to obtain 

Hy|li 2 2(r) +^P HVr(t)2/||| 2(r(t)) + J ll2/ll^(r W )dt<c(||y||^ 1(r(0)) + ||/||^ 2(r J . (3.10) 

[0,T] 

Note that from [LM68, Ch. I,Thrm. 3.1] it then follows that <fty G C([0, T], i/ 1 (r(s))). 
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4 Control constrained optimal control problems 

Using the results from the previous section, we can now formulate all kinds of control- 
constrained optimal control problems known for stationary domains, see for example [TV605]. 
As a first example, given a moving surface as in Assumption 2.1, let St '■ L 2 L2 ^ — > L 2 (T(T)) 
denote the solution operator u i-> y(T), where y satisfies 



J yfdT(t)+ J V r y • Vr^dr(t) = (0,y) H -i (rW)iH i (rW) + (n,v3) L 2 2 ^ , (4.1) 



dt 

r(t) r(t) 



for all ip G Wr, and with y(0) = G L 2 (ro). We know, that every function y G Wp has a 
representation in C([0, T], L 2 (T(s))) for any s G [0, T], compare Lemma 2.2, and the inclusion 
4>(.jWr C C([0, T], L 2 (r(s))) is continuous (in fact compact). Thus St is a continuous linear 
operator. Consider the Control problem 



f min ueL 2 2(r) J(u) := |||5 T («) - 2/r || | 2(r(T)) + f IML^ 
s.t. a < u < b , 



with a, a, 6 G M, a < 6 , a > 0, and yT G L 2 (T(T)). This is now a well posed problem. By 
standard arguments, see for example [Tr605, Thm. 3.15], using the weak lower semicontinuity 
of </(•), one can conclude the existence of a unique solution u G L 2 L2( j.y 

For an other example let the linear continuous solution operator S d : L 2 L2 ^ — > L 2 L2 ^ r y u i-> y, 
where y solves (4.1), and consider the problem 

i 2 (r) L 2 (r) i2 ( r ) 

S.t. a < U < b , 

with q, a, b as above and G L 2 L2( j,y Again there exists a unique solution, see [Tr605, Thm. 
3.16]. 

The first order necessary optimality condition for (P d ) reads 

(S d u-y d ,S d (v-u)} L 2 +a(u,v-u) L 2 = (au + S* d {S d u- y d ),v - u) L 2 > , (4.2) 

L 2 (r) £ 2 (r) L 2 (r) 

for all v G = G L 2 h2 ^ \ a < v < 6 j. The adjoint operator S*^ : L 2 h2 ^ — >■ L 2 L2 ^ maps 
w G L 2 L2 ^ onto the solution p G Wr of 

-(f 1 , ¥») j?-i(r(t)),ifi(r(t)) + y V r p • V r ^ dr(t) = (v, ^) £ a ^ , (4.3) 

r(t) 

for all 99 G W r , and p(T) = G L 2 (r(T)). This follows if one tests (4.1) with p and (4.3) with 
y. Integrate over [0,T] and use y(0) = and p(T) = to arrive at (v,y) L 2 = (p,u) L 2 , 

1,2 (r) L 2 (r) 

for u,v £ ^^2(- r -j arbitrary. 

Note that via the time transform t' = T — t Equation (4.3) converts into equation (3.8) with 
b = — divp^F. Therefore all the results from Section 3 also apply to (4.3). 
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The necessary condition (4.2) characterizes the optimum u as the orthogonal projection of 
~h^d(^d,u — yd) onto Uad. In our situation this is the pointwise application of the projection 
P [aM :R^[a,b]. 

Lemma 4.1. Let Pu ad denote the L 2 L2 ^ T y orthogonal projection onto U a( i, which is defined by 
{u-P UaA {u),v-P UiiA { U )) Ll2 ^<U VveU ad . (4.4) 

Then for u G L 2 L2 ^ one has for a.e. t G [0, T] 

Pu ad (u)[t] =P M («[*]). 

Proof. Because L 2 ([0, T], L 2 (r(s))) is isometrically isomorph to L 2 ([0, T] x T(s)), and because 
<t>fP [aM (u) G L 2 ([0,T] x 1»), we also have <f> s t P [aM (u) G L 2 ([0,T], L 2 (T(s))). ' 
Let C = {t G [0,T] | Pf/ ad (u)[t] / J P[ a ,fe](n[t]) } and assume meas(C) > 0. Now test (4.4) with 
v = P[ a ^(u[t]) to arrive at 

/ («-^aaM.^]W-^W>^(r (t ))di<0. (4.5) 

[0,T] 

But now for a.e. t G [0, T] and a.e. 7 G T(t) one has 

KM - pM^MX^mMH - % a >t)[7] > 0, 

because -P;/ ad (wt)[7] G [a, 6]. Moreover for t G C we have 

KM - p^OMX^mWH - %ad(^)H > 0, 

on a set of positive measure. Since meas(C) > this contradicts (4.5). I 
Introducing the adjoint state Pd{u) = S d (SdU — yd), let us now rewrite (4.2) as 

u = P [aM (-^p d (u)Y (4.6) 

Similarly the unique solution u of (Ft) is characterized by u = P[ a ,b] {—^Pt(u)) , with pt(u) = 
S t (Stu — |/t)- Note that however the adjoint state pt in general is less smooth than pd- This 
is because the adjoint equation, i.e. the equation describing S T : L 2 (r(T)) — > L 2 2( - r ^, v ^ p, 
reads 

-{p,f)H-^r(t)),m(r(t))+ J V r p • V r ^dr(t) = , 

r(t) 

for all if G Wr and with p(T) = v G L 2 (r(T)). While Theorem 3.10 applies, this is not the 
case for the smoothness assertion (3.10), as long as yd G L 2 (r(T)) \ P 1 (r(T)). 
Before we can discuss the discretized control problems in Section 7, in the next two sections 
we present some results on the discretization of the state equation. 
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5 Finite Element Discretization 



We now discretize T using an approximation Tq of To which is globally of class C 0,1 . For the 
sake of convenience let us assume n = 2, i.e. T(t) is a hypersurface in M 3 . 
Following [Dzi88] and [DE07], we consider Tq = \J ieI consisting of triangles T % h with 
corners on To, whose maximum diameter is denoted by h. With FEM error bounds in mind 
we assume the family of triangulations {Fq}/j >0 to be regular in the usual sense that the 
angles of all triangles are bounded away from zero uniformly in h. 

As detailed in [DE10] and [DE07] an evolving triangulation T h (t) of T(t) is obtained by 
subjecting the vertices of Tq to the flow <!>. Hence, the nodes of T h (t) reside on F(t) for all 
times t G [0, T], the triangles T^ being deformed into triangles T^(t) by the movement of the 
vertices. Let rrih denote the number of vertices {Xj}™ 1 ^ in Tq. Now Xj(t) solves 



d 



dt X j (t) = V(X j (t),t), X j (0) = X°. (5.1) 
Consider the finite element space 

Vfr(t) = [v G L 2 (T h (t)) | v G C(T h (t)) and Vi G I h : G U l (Tl(t)) } 

of piecewise linear, globally continuous functions on T h (t), and its nodal basis functions 
{(^•(i)}™^ that are one at exactly one vertex Xi(t) of T h (t) and zero at all others. For the 
finite element approach, it is crucial for the triangles T£(i) not to degenerate while T h (t) 
evolves, which leads us to the following assumption. 

Assumption 5.1. The angles of the triangles T^it) are bounded away from zero, uniformly 
w.r.t. h, i and t. Also assume at(T h (t)) = T(t), with the restriction of a t to T h (t) being a 
homeomorphism between T h (t) and T(t). 

In order to ensure optimal approximation properties of the discretization of the surface, we 
require d to be twice Lipschitz-continuously differentiable. 

Assumption 5.2. d G C 2 > l {N T )- 

Let us summarize some basic properties of the family {^ h (t)}te[o,T]- 

Lemma and Definition 5.3. Let <& s h : T h (s) x [0,T] ->■ R 3 denote the flow ofT h , i.e. the 
unique continuous map, such that $| h (7^(s)) = T^t) and <& s th is affine linear on each T^(s). 
There holds <& r t h = <& s t h o $ r h and thus & s h o ® s t h = id T h( s y The velocity Vh = d t ^ h is the 
piecewise linear interpolant ofV on each triangle T^it). 

As in Lemma 3.2 we define the pull-back 4> s th : L 2 (T h (t)) ->■ L 2 (T h (s)), <j) s th v = vo& sh . 
The piecewise constant Jacobian J* h of $>f h satisfies for all s,t G [0, T] 

-4 < min J- h ( 7 ) < max J« ft ( 7 ) < C h 3 , (5.2) 

for some constant Cj > that does not depend on h > 0. Moreover J* h and Dp( s )$f ft : 
TF(s) — > TT(t) C M 3 are differentiable with respect to time in the interior of each T^{s). 
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The nodal basis functions have the transport property 

Vi = 4>o,h^ < f>t,h ( Pi = > ^<i<m h . (5.3) 
Let u h (t) denote the normals ofT h (t), defined on each T^(t). 

Proof. Consider a Triangle T l h (s), s G [0,T]. W.l.o.g. let Xi(s), X 2 (s), X 3 (s) denote its ver- 
tices. Then, using matrices X l (t) = {X 2 (t) - Xi(t),X 3 (t) - Xi(t)), we can write 7 G T£(s) 
in reduced barycentric coordinates as A 7 (s) = (X l (s) T JP(s)) _1 X ? (s) T (7 — Xi(s)). On T£(s) 
the transformation is uniquely defined by A$» h7 (i) = A 7 (s) and thus 

^(7) = X i (t)(X i ( S fX i ( S ))- 1 X i ( S ) T ( 7 - Jfi( a )) + Xi(t) . 

In the relative interior of T^{s) the map & s th '■ THs) —> T^(t) is differentiable and its derivative 
T) T i^<& s th : R 3 d TT^(s) — > TT^(t) C IR 3 can be represented in terms of the standard basis 
of M 3 by the matrix Dj jt = X i (t)X i (s) T X i (s))- 1 X i (s) T . 

Now one easily proves that the angle condition in Assumption 5.1 ensures the existence of c > 
such that \ T X i (s) T X i (s)X > c min(||X 2 (s)-Xi(s)|| 2 , ||X 3 (s)-Xi(s)|| 2 )||A|| 2 for all A G R 2 , 
s G [0,T]. Hence, \\(X*(s) T X* (*))"% < (c min(||X 2 ( S ) - Xi( S )|| 2 , \\X 3 (s) - X^f))' 1 , 
and since ||X'(s) T ||l, ||X*(s)||| < 2max(||A 2 (s) - Ai(s)|| 2 , \\X 3 (s) - Ai(s)|| 2 ) we get 

\\ B n(s)*!, h ^\\ < c _max(||X 2 ( S ) - X 1 (s)\\ 2 , \\X 3 (s) - X 1 (s)f) 



||d 7 || " min(||X 2 ( S ) - X 1 (s)\\*, \\X 3 (s) - X 1 (sW 

for all d7 G TT^(s). Using again Assumption 5.1 one concludes that the quotient of edge 

lengths is uniformly bounded. 

Also, one easily verifies for r,t G [0, T] 

$t,h7 = (*t,h C,/J7 and $^$ t % = id rh(s) . (5.4) 



y det ((D T i (s) $* J T D T ; (s) $*J on the triangle T l h (s), where the deriva- 
tive is represented with respect to an orthonormal basis 03 (s) of TTHs). As per above consid- 
erations the spectral radius of {D\ t ) T D\ t is uniformly bounded. Hence, there exists Cj > 



We have J^ h 



such that Jf h = < /det (*B(s) T (.D* t ) T ^ > l t^( s )) — Cj. Because we can switch s and t 
and since by (5.4) we have ((fr^) -1 = <& l sh and thus ji- = J l s h < Cj we conclude 

Vs ,tG[o,T] : v 7 erji< Jl h in)<c h j. 

The trajectories $^7, 7 G T h (s), the Jacobians J^, and the entries of D l s t are differentiable 
for t, because the trajectories Xj(t), 1 < j < are, compare (5.1). Hence also D r ^ s ^f h 
is differentiable as a map into R 3 . The velocity ^(7, s) = dt$ s t fe 7 equals V at the vertices 
and depends linearly on the coordinates A 7 . As for the transport property (5.3), it is a 
consequence of the piecewise linear transformations of the piecewise linear Ansatz functions 
ifi which implies $ h (fi{t) = (fi(0), compare [DE07, Prop. 5.4]. I 
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Remark 5.4. Similarly one can prove the map <& s th : T h (s) — > T h (t) to be bi-Lipschitz with 
respect to the respective metrics. The Lipschitz constant L does not depend on s,t G [0,T]. 

In order to compare functions defined on F h (t) with functions on T(t), for sufficiently small 
h > we use the projection a t from (2.1) to lift a function y G L 2 (T h (t)) to T(t) 

y\a t (x)) = y{x) Vx G r ft (t) , 

and for y G L 2 (r(t)) we define the inverse lift 

m (x) = y(a t (x)) Vx G T h (t) . 

For small mesh parameters h the lift operation (•); : L 2 (T(t)) — > L 2 (T h ) defines a linear 
homeomorphism with inverse (•)'. Moreover, there exists Cmt > such that 

1 - c illt h 2 < ||(-)Hl£(L 2 (r(t)),L 2 (r h (t)))' IK')'ll£(L 2 (r h (t)),L 2 (r(t))) - 1 + c mt^ 2 , (5-5) 
as shows the following lemma. 

Lemma and Definition 5.5. The restriction of at to F h (t) is a piecewise diffeomorphism. Denote 
by 5 h the Jacobian of a t \ V h {t) : T h (t) T(t), i.e. 6 h = ^ = |det(M)| where M G M 2x2 repre- 
sents the Derivative da t (x) : T x T h (t) —t T a ^T(t) with respect to arbitrary orthonormal bases 
of the respective tangential space. For small h > there holds 

sup sup |1 — 8h\ < Ch 2 , 

te[o,T] r(t) 

In particular at\i>h( t } is a diffeomorphism on each triangle T^t) . Now ^- = = |det(M _1 )|, 
so that by the change of variable formula 

< Ci llt h 2 \\v || L i (r) . 

Also there exists C > such that 

1. sup t6 [ 0j T] II^C^) ll_L oc (r^(t)) ^ Ch 2 , where the material derivative is to be understood in 
the sense of h and 

2. sup te[0jT] \\V{I - n l h )V\\ L °o {m) < Ch 2 , where K h = ^ (/ - dU) V h (I - dU), U i3 = 

d Xi Xjd, and V = {5ij — viv 'jYtj=\ an d 7~ >h = {<% ~ u i u j}7t=i are P ro i ec ^ ons on the 
respective tangential space. 

Proof. We summarize the proof given in [DE07, Lemma 5.1] to extend it for the 1. assertion. A 
similar proof can be found in [DE10, Lemma 5.4]. Following [DE07], we use local coordinates 
on a triangle e := T^(s) . W.l.o.g. one can assume e G M. 2 x {0}. Since both d and d = ji<f>t 
equal zero at the corners, the linear interpolates Ihd, l^d vanish on e thus, using standard 
finite element approximation results, we get 

NIU<»(e) = \\d~ hd\\ L ^(e) < C/J 2 ||d||/P,°°(e) < ch 2 \\ d\\ C l,l {Nt) 



^dr ft (t) - J vAT{t) 
!•'■(/! r(t) 



61 



vdT(t) 



r(t) 
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and similarly ||d||L°°(e) < c ^ 2 ||^IIc 2 > 1 (.a/'t)- Also one has 



\dxid\\ L oo {e) < c/t||d|| c i,i(A/- T ) and ||a Xi d||£oo( e ) < ch\\d\\ C 2,i (A f T) 



(5.6) 



for i = 1, 2 at any point (xi, x 2 , 0) G e. 

Consider the basis Q5(i) = {c^^ h, d X2 ®thi ^(t)} of M 3 , whose first two members span the 
tangential space of T^(t). Let (y\(t), 1/2 (t) , u^t)) represent = Vd(-,t) with respect to 
*B(t). Note that 05 (s) are the unit vectors. 

We have (vi(t), z^ 2 (t)) T = M t -1 (D( Xl X2 ^^ h ) T Vd, with the uniformly positive definite matrix 
M t = (B {xiiX2) ^l h fD {xuX2) ^ h . Now 

0{h) = D {xuX2) d' = B (xijXa) ^</> a tth d =^>{ Xl , X2 )<t> s t , h d = ^((j)l h Vd T D (xuX2) <S>l h 



d/ 



d/ 



©(h) 



where we used ^(7, s) = d Xt d(^, s) and (5.6). We subsume 

INIU 00 ^); INIU 00 ^) < c ^IMIIc* 2 . 1 (A't) • 

One has 

Dot = Id - Vd(Vd) T - dV 2 d 
and with Vd(-,s) = (^i(s), f 2 (s), ^3(s)) T we compute (see [DE07]) 



^ = ||^iatx5x2«t|l = l^l+d-R^dzi^dra^) = V 1 - ^1 - vl+dR{v,d Xl v,d X2 v) = l+0(h 
with some smooth remainder function i? and finally, since <i = 0(h 2 ) 

-V\i>\- f 2 Z> 2 



II^IU°°(e) 



^1 - ^ 2 - ^ 



+ C(/i 2 



L°°(e) 



where we used |z>j| < CTi, i = 1,2. For a proof of 2. see [DE07, Lemma 5.1]. I 
The next Lemma concerns the continuity of the lift operations between L 2 L2 ^ h ^ and L 2 L2 ^ r y 

Lemma 5.6. Using the pull-back 4>\ h we can define L 2 2 ^ r ^ as in Lemma 3.4- For sufficiently 
small h > the lift operation (•)' constitutes a continuous isomorphism between L 2 L2 ^ and 
L\ 2 (T h ) Averse (■)/. There holds 



L 2 {V h ) l 2 (t) 



< c iQ th 2 \{f,g) L 2 



L 2 (T)' 



Proof. Let L 2 L2 ^ Ti y according to the flow ^ h as defined in Lemma 3.4. We define L 2 



J L 2 {V h ) 



{Jiel h ^ L 2pn) with the scalar product f Q (-, ■) L 2^ T h^dt. 

Let ft = ^QOOfOiJ)^ denote the mapping between Tq and To induced by the projection at- 
By Assumption 5.2 and by the construction of <J>£ and is follows that f t : Tq — > Tq is 
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a diffeomorphism on each triangle T^(O) and globally one-to-one and onto. Also ^ and its 
spatial derivatives are continuous w.r.t. time t. 

We will show that I : Tq x [0,T] ->• T x [0,T], (7,*) (* t ( 7 ),f) is a piecewise diffeo- 
morphism whose Jacobian is bounded away from zero. By Assumption 5.1 we already have 
that * is globally one-to-one. Together this implies that the pull-back with ^ constitutes an 
isomorphism between L 2 (T Q x [0,T]) and L 2 (Tq x [0,T]). This again means that 

G L 2 ([0,T],L 2 (rg)) ^ G L 2 ([0,T],L 2 (r )) . 

As to * being al local diffeomorphism, the sets T l h = \J te [ T j T£ (i) are a partition of Tq x [0, T] . 
In the interior of each the map 'I' is a diffeomorphism. In fact, let 7 G int(T^) for some 
1 < z < m/j. Compute 



D 



rgx[o,T] 



*(7) 



D r& * t ( 7 ) ^(7) 
1 



We have Dph^j = Dp^^gDph^-jafDph^^. Its Jacobian is the product of the Jacobians Jg, 
5^, and J® h that are each bounded away from zero, uniformly in 7 and t, compare (5.2), and 
the Lemmas 5.5 and 3.2. Hence the Jacobian of ^ is bounded away from zero. 
As to continuity of (■)/, by Lemma 5.5 we have that 



(fl,9i) L 2 9 , - (/ ) 5 , )l2 



0,T] r(t) 



fg{^-l)dT{t)dt 
°h 



< Cinth 2 \{f,g) L 2\ 



Now, instead of dealing with Problem (3.8) directly, w.l.o.g. we consider the equation 

^ J yipdT{t) + J V r yVr^p + tJ-y(pdT(t) = {ip 1 y) L 2 (m) + {f,ip) L 2 {m) , (5.7) 



r(t) 



r(t) 



with p, G M. large enough to ensure ji := b + p, > 1. Note that y solves (5.7) iff e^y solves 
(3.8) with right-hand side e^f. 

In order to formulate the space-discretization of (5.7), consider the trial space 



m h 



H k = £ j/i(*)¥>i(*) e L\ 2{Th) I G ^([0,T]) ~ H\[Q,T\) m » . 



The following definition of weak material derivatives for functions in Hy exploits the fact 
that Hy is isomorph to H 1 ([0,T]) mh . It thus avoids the issue of extending the theory from 
Section 3 for the smooth surfaces T(t) to our Lipschitz approximations T h (t). 

Lemma and Definition 5.7. The weak material derivative of v = Yli^iViifytPiify e ^v h 
is v = (pQ h {<j^l h v)' = Yl^i v'i{'t)<Pi{i) ■ Let further w G Hy , then (v, w) L 2^p h (t)) * s absolutely 
continuous and 



d_ 
dt 



J vwdT h (t)= J vw + vw + vw dwr h V h dT h (t) . 



i*(t) 



i*(t) 
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Proof. Observe v = <^(<^>)' = <j>* th (ETA «i(%i(0))' = (££\ 3(%i(0)) because 
(^ )h v(t))'(7) = 5W(0)(7) = for all 7 G I* as in (5.3). 

Apply Lemma 3.1 on each triangle to see that (<Pi(t), <Pj(t))L 2 (r h (t)) 1S smooth and 

-^{<Pi{t),(pj(t)) L 2 {vh{t)) = I ipwj dW Vh V h dT h (t). 

i*(t) 

Now 

{v,w) L 2 {rh{t)) = ^ ^(*)^(*)(^(*)^j(i))L2 (r / J ( t )) 

and the second assertion follows, since Vi,Wj G i7 1 ([0, T]), 1 < i,j < m^. I 

We approximate (5.7) by the following semi-discrete Problem. Consider a piecewise smooth, 
globally Lipschitz approximation A of m, such that A > 1. Find y G Hy h such that for all 

^ J y h LpdT h (t)+ J V r hy h -V r h(p+Xy h <pdF h (t) = {ip,y) L 2 {r h W) + {f h ,<p) L 2 {Th{t)) , (5.8) 
r h (t) r h (t) 

and y h (0) = yft G 14(0). One possible choice would be A = fj, h f h = ft and yfi = Pft((y )i) 
with Pq the L 2 (rQ)-orthogonal projection onto Vh(0). 

First of all let us state that (5.8) admits a unique solution in H v . This is because for 
yh = Yli=i Vi^Pi we can rewrite (5.8) as a smooth linear ODE with non-smooth inhomogeneity 
for the coefficient vector y = {yi}^ G fl^QO, T]) mh 

±(M(t)y h (t)) + (A x (t))y(t) = F(t), y h (0) = y£ , (5.9) 

with smooth mass and stiffness matrices 

M(t) = {{w, ^> L 2 (rftW )}™/ =1 and A\(t) = { J V r H<PiV r H<Pj + \<Pi<Pj dT^*)}™^ , 

i*(t) 

and right-hand side F(t) = {(fuVi) tf(Th{t))}7=i G L 2 ([0, T], R mh ), compare also [DE07]. 
Observe that we used the continuity of the coefficients £/j G -ff 1 ([0,T]) as well as tpi = 0. 
Existence of a solution G -ff 1 ([0, T]) mh of (5.9) can be argued by variation of constants or, 
more generally, one can apply an existence result by Caratheodory, compare [CL55, Thms. 
1.1+1.3]. Uniqueness of yh is a consequence of the following lemma. 

Lemma 5.8 (Stability). Let yo G L 2 (Fo) and f G L 2 L2 ^ v y and let yh solve (5.8) with y$ G 
Vh(0) and fh = fi- There exists C > 0, such that for sufficiently small h > the solution 
satisfies 

\\Vh\\h(rk(T)) + J J ^r^h) 2 + W h dT h {t)dt < C (\\y%\\ 
as well as 



2 

L h ( T) 



\y h f +ess sup f (V Th y h f + \yldT\t)<c(\\y%\\ 2 H , {T H ^WfWh 
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Proof. From the definition of M and A\ using Lemma 3.1 on each triangle T^(i) there follows 
M'(t) = { J ^div r ^dr ft (i)}™/ =1 and 
r*(t) 

d. /* 

(^40y = / -Vrfc^i (D rh T4 + D r hVf ) V rh (pj + A^^j + . . . 
r*(t) 

+(V r h(/3jV r h(/9 : ,- + \ipi<pj) div F h 14 dr ft (t) . 

Multiply (5.9) by y' to obtain 

ll*"»ll£ 2(r h (t)) 

y W + ^ (yA A y) = -y'M'y + ^yA' A y + Fy' 

< c( lly?illffi(r(t)) +H/'llL2(r h (t))) + 9 II ^ II L^r^ft)) ' 

\ v v s ' ^ 

and a Gronwall argument yields the second estimate. Multiply (5.9) by y and proceed similarly 
to prove the first. I 

Obviously the material derivative depends on the evolution of the surface, i.e. different deriva- 
tives arise according to whether 4>f or c/>| h is applied to pull back a function to a fixed domain. 

In order to compare z l h with (zhf we need the following lemma. 

Lemma 5.9. Let y = Y^i=i Vi^Pi £ ^v h - The lift y l lies in Wr with y l G L 2 L2 ^ r y and for a.e. 
t G [0, T] there holds 

y l -(y) l \ <Ch 2 \\V m y l \\, 

a.e. on T(t). 

Proof. We start by computing the material derivatives of (f>i(x,t) : Nt — > <f>i(x,t) = 
if\(at(x),t), i.e. the constant extension of the trial function ipi, 1 < i < rrih, along the normal 
field of r(£), compare the proof of [DE07, Thm. 6.2]. Observe that tp\ is not smooth along the 
edges of patches at(T^(t)). However, <p\ is smooth in the (relative) interior of all at{T ] h (t)) . 
Derive (pi at 7 € remit (TjJ(i)) to obtain 

V^( 7 , t) = V^(a t ( 7 ), t) (Id - Vd( 7 , i)Vd( 7 , t) T - din, t)V 2 d( 7 , t)) 

^ i (7,t) = ^(a t (7),t) + V^(a t (7),t)(-^d(7,t)Vd(7,t)-d(7,t)5tVd(7,t)). 

By construction of (pi we have Vy?j(aj(7))Vd(7, t) = Vr</?'(ai(7)) W(at(7), i) = since (pi is 
constant along orthogonal lines through T. Also, from d(<&° (7), i) = it follows <9td = — VdV. 
The (strong) material derivatives do not depend on the extension (p~i, but only on the values 
on T and T h , respectively. One gets (p\(a t (i),t) = d t <fi(a t (-f),t) + VCpi(a t {^),t)V(a t (~f),t) and 
ipi(~/,t) = dt(fi(^,t) + V (p~i(n ,t)Vh(n ,t) which together with (5.10) leads us to 

v\ = +(V-V h + d((V 2 d)V h + d t Vd)) V m <p\ , (5.11) 
in the relative interior of the patches at(T^(t)), j € 1^. 
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In order to prove that the pull-back p := $p\ lies in C 1 ([0, T], L 2 (r )) D C([0, T], H 1 ^)) 
for all 1 < i < rrih we proceed in four steps. 

1. We show that p is globally Lipschitz on Tq x [0,T]. Observe, that (5.10) implies that all 
derivatives of p exist and are bounded on the interior of patches -P^(i) = $g(at(Tj£ (t))). Since 
$i = $5 o fl t o : rj x [0,T] ->■ T smoothly maps the edges of Tq into T the domains 
Ute[oT] x W C To x [0,T] have piecewise C 1 -boundaries. Also, </5 is continuous and we 
are in the situation to apply Stoke's theorem to confirm p G W 1,00 (To x [0,T]). By Morrey's 
lemma, for a formulation on manifolds see [MMV98], we conclude <p G C 0,1 (ro x [0,T]). 

2. Now as to the time derivative, fix e > and t G (0, T). Let L > denote the global Lipschitz 
constant of p on To x [0, T] and choose 77 > sufficiently small such that Yliei h meas(P^ \ 
Pl. q ) < e 2 /8L 2 where P l hr) = {7 G P£ | ^(7) C P^}, the balls ^(7) being taken with 
respect to the metric of To- Now, as stated above, the patches Pl(t) = &(t)(Tl) move con- 
tinuously across To, and we can choose K sufficiently small such that for all i G Ih and 
k G (-K, K) we have P^ v (t) C P l h (t + k). The derivative d t p(^,t) = $(p\ which is defined 

a.e. on To x [0, T] then is continuous on the compact set )C V = Uie/ h ??(*) x [t — K,t + K] 
and we have 

p J m + k) - (p{t) - d t p{t)k) 2 dr = p ( / + fe ) - - 5 ^W fc ) 2 dr o 



+ J {<p(t + fc) - <p{t) - d t p{t)k) 2 dr 

Substituting ^(7, i + fc) — ^(7, i) = d t (p(j, t)k + J^{d t p{^, t + rk) - d t p(j, t))kdr on P l h t like 
in the proof of Lemma 3.6 we choose k small enough for 

sup \\dMt + rk)-d t p(t)\\l< "l (v v (5-12) 
re[o,i] 2meas(l J 

which is possible by uniform continuity of dtp on 1C V . Estimating the second addend by 
{ZLk? J2 ieIh meas(P* \ P^) < e 2 /2 yields 

limsup-||£(* + k) - <p(t) - d t <p{t)k\\ L 2 (ro) < e. 
fc^O K 

for every e > 0. Hence p is differentiable into L 2 (Fq) with derivative dtp. 
3. Thus in order to show p G C^QO, T], L 2 (T )) it remains to prove that d t p : [0,7] -> L 2 (r ) 
is continuous. By (5.10) is essentially bounded on To x [0, T]. Let M = ||<9tV?||,L°°(r x[o,T])- 
For e > choose rj > sufficiently small such that J2i£i h meas(P^ \Ph v ) — e 2 /8M 2 . As 
above, choose if > and /C,, accordingly. Now, choosing A; > small enough such that (5.12) 
holds one arrives at 

\\d t p(t + k)-dtp(t)\\h { r ) = T,( I (dtp{t + k)-d t p(t)) 2 dT ... 



+ J (d t p{t + k) - d t p{t)) 2 dT^j < e 2 
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4- Continuity of p : [0, T] — > H 1 (Tq) follows similarly. In fact, the spatial partial derivatives 
of <p exhibit the same piecewise smooth structure as dtp. 

Finally, (p = $<p\ G C^flO, T], L 2 (F )) n C{[0,T],H r ) implies G W o (0,T), and we 

conclude y l G Wr as well as y G ^2/ r y The estimate now is a consequence of (5.11). I 

Before we proceed to the main result of this section, we need to understand the approximation 
of elliptic equations on T(t) by finite elements on T h (t). 

Lemma 5.10. For t G [0,T] and g G L 2 {T(t)), g h G L 2 (T h {t)) consider 

J v r z° ■v r p + fiZVpdr(t) = (g,p) L 2 im) ,VpeH 1 (r(t)) (5.13) 

r(t) 

and 

J V r „zf ■V Th cp + fi l Z 9 h h l pdT h (t) = (g h ,i P ) L 2 {rh{t)) ,y<peV h (t) (5.14) 
r*(t) 

with unique solutions Z 9 G // 1 (r(t)) and Z^ G V^(t). TTie solution operators S(t) : L 2 (r(i)) — > 
L 2 (r(t)), 5 .-»• Z» and S fc (t) : L 2 (r fe (t)) VJ, C L 2 (r fc (t)), g h ^ ZjJ h are self-adjoint. There 
exists C independent of t G [0, T] suc/i that 

1. \/p G Vfe(t) : 1 \\ip l ||^i(r( t )) - IM^r^t))! ^ Ch2 \\ l P l \\m(r(t)) < 00 as weU as 

2. \\(-) l S h {t)(-f - S(t)\\ c{LHmhLHm)) < Ch 2 and 

3. \\(-) l S h {t)(-T - S(t)\\ c{LHrmm(rm < Ch. 

Proof. The operators being well-defined and self-adjoint follows by standard arbuments. As- 
sertion 1. follows from Lemma 5.5 [2.], since ip l is continuous and piecewise smooth on T(t) 
and thus lies in i? 1 (r(t)) with 

V rh p\\ 2 dT h (t) = J \\V r p l \\ 2 dT(t)+ J VrV (n l h - Id) VrV dT(t) , 
r*(t) r(t) r(t) 

for details see for example [DE07, Lem. 5.2] and proof. 

For a proof of 2. and 5. see [Dzi88, Thm. 8] and the discussion of (•); and (-) z preceding 
Lemma 4 in aforementioned article. The fact that C does not depend on t is a consequence 
of Assumption 2.1 and 5.1. I 

Theorem 5.11. Let Assumption 2.1, 5.1 and 5.2 hold and let y G Wr solve (5.7) for some 
f G L 2 L2( j,y yo G ^ 1 (ro), such that (3.10) holds. Let y^ solve (5.8) with A = pi and fh = fi 
and some approximation y$ of (yo)i- There exists C > independent of y and h such that 

\\v l h - y\\l lHr) < c (\\ Vh (o) - m(o)\\l Hrho) + ^ (llyo||^ (ro) + + H/%J) • 
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Proof. Define z = S(t) (y l h — y) and = Sh(t) (5h (yh — Vi)) with S(t) and Shit) as in Lemma 
5.10. Now 5h (yh — Di) = (•)'* (4 ~ v) an ^ hence it follows from Lemma 5.10[#.] that 

114 - 4m m) = ll(0)^(-f - S)(y{ - y)\\ LHm) < Ch 2 \\y l h - y\\ L 2 {m) , (5.15) 

Observe now for Zh = Yli^i Zifi using Lemma 5.9 we get 

Y = {(y[ - y,^ L 2 {m) }Z\ e H\[0,T\) m » , and thus z = (A^Y G H\%T])^ . 

Hence z G Hy h and again by Lemma 5.9 z l h G Wr as well as z l h (t) G L 2 (T(t)). 
We can now test (5.7) with z l h , using (5.13) in the process, to obtain 

^ {y, z l h)mv{t)) + (y, 4 - y)^(r(t)) = (4> ^(rw) + </> 4>L 2 (r(t)) • • • ^ 5 16 ^ 

+ {-A r y + fiy,z- z l h ) L 2 im) , 

and testing (5.8) with z^ gives 

^j.(yh, z h)L*(rh(t)) + (4> 4 - y)^{v(t)) = {zh, yh)^^^) + {fu *h) . (5.17) 

Now, since the strong material derivative exists and is continuous on each triangle T^(t), 
the scalar products (ipi, <Pj&h) L 2 (T h (t))i 1 — ' — m ^> are differentiable with 

i*(t) 

and we have 

~j^4> 4)i 2 (r(t)) =-^(yh, z hbh)L 2 (T h (t)) 

= dt^ yhl Zh ^ L2 ( rh ( t )) + ^ ~~ 1 ))i 2 (r' i (t)) + (y/n ^^)L 2 (r h (t)) • • • 
+ (yh,z h (Sh - l)>L 2 (r h (t)) + (Vh,Zh div T hV h (S h - l)) L 2 {r h (t)) . 

Hence, we can rewrite (5.17) by means of the L 2 (T(t)) 

^(vl 4>L 2 (r(t)) + (yl y{ ~ y)^(r(t)) = ((4)', 4>L 2 (r(t)) + if, 4>L 2 (r(t)) + R h , (5.18) 
with 

R h ={yh,z h 5 h ) L 2 {rh{t)) + {yh,z h (Sh - l))L 2 (rh(t)) + (Vh,Zh div rh V h (5 h - l)) L 2 (r h W) ... 
+ (fi, zh{l - $h))L 2 (rh(t)) ■ 

Subtracting (5.16) from (5.18) yields 

^(4 ~ V, 4>L 2 (r(t)) + \\y l h ~ y|li 2 (r(t)) = ((4)' - 4> ^ 2 (r(t)) + ^h, (yh ~ yi)5h) L *{rh{t)) ■ ■ ■ 

+R h + (-A T y + py, z l h - z) L 2 (m) . 
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From (5.14) we know (z h ,S h {y h - yi)) L 2 {V h {t)) = z' h A x z h = \^ t (z h A x z h ) - \z h A' x {t)z h , in the 
notation of (5.9). Now, using (5.15) and 



\R h \ < C7i 2 ||z fe || L 2 (r fc W ) (||y/ l ||x,2( r h( t )) + \\yh\\L 2 (r h (t)) + II/jIIl 2 ^' 



(*)) 



we can estimate 



Id 

2dt 



(z h A x (t)z h ) + \\y l h - y|||2( r(t )) <C^ 2 ||y||L 2 (r(t))l|V r h(t)Z/ l ||( L 2( r h W ))n=i . . 



+ \\ z h\\ 2 Hi(r h (t)) + h 2 \\y\\H 2 (r(t))\\yh - ^11^(^(4))) + \R h \ 

<^hh - yiWh^hit)) + c(z h A x (t)z h . . . 

+ /l4 (felli 2 (r'>(t)) + ll&"illz, 2 (rfc(t)) + Wfi\\ 2 L 2 (r h (t)) + IMI# 2 (r(t)))^ • 
We can now apply Gronwall's lemma for 

[z h A x {t)z h ]l + J ||j/ft-j/||i 2(r(t)) d* 

< C7i / ll?//i|li,2(rh(t)) + lly^llL 2 (r h (t)) + II /ill £2(^(4)) + lbllif 2 (r(t)) , 

[0,T] 

and with the stability estimate (3.10) and the Lemmas 5.8 and 5.6 we finally arrive at 

=(y l h (o)-y(o),z h ) L 2 (ro) 



J (v rho z h ) 2 + \z 2 h dT h ... 



(5.20) 



+ h 4 ^||z/o|lii ( r ) + ho H^i(rg) + 



L 2 (r) 



Apply again (5.15) to prove the lemma. 



Remark 5.12. Depending on the regularity of yo, possible choices of yft yielding 0(h?)- 
convergence of y l h comprehend the piecewise interpolation of (yo); an d the L 2 (ro)-orthogonal 
projection of (yo)/ onto Vh{0). For the latter, the term involving Zh in (5.20) vanishes com- 
pletely, but it's i? 1 (ro)-stability requires further investigation. 

The order of convergence is lower, if the solution of (5.7) does not satisfy the additional 
regularity estimate (3.10). 

Theorem 5.13. Let Assumption 2.1, 5.1 and 5.2 hold and let y 6 Wr solve (5.7) for f = 0, 
and y G L 2 (T ). Th ere exists C > independent of y and h such that for the solution y^ of 
(5.8) with yo = Po({yo)i) and fh = there holds 

\\y l h -y\\ 2 L 2 <C[h 2 + sup || A' -Hli~(r(t)) J IIJA)llia(r ) ■ 
i 2 (n V te[o,T] J 
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Proof. We proceed as in the proof of Theorem 5.11 up to (5.16) which now reads 

^(y, 4>L 2 (r(t)) + {y, y[ ~ v)tf(r(t)) = (4> ^ 2 (r(t)) + (~ A ry + W,z- z l h )H-HY{t)),m(Y(t)) , 

Analogously to (5.15) we can apply Lemma 5.10[3.] and estimate the last term through 

\{-A r y + Xy,z l h - z) H -i {m)jH i {r{t)) \ < \\ - A r y + Aj/|| H -i( r(t ))||4 - 4m(Y(t)) ■ ■ ■ 

< || - A r y + Aj/|| H -i( r(t ))C/i||^ - y\\ L *(T{t)) ■ 

On the other hand (5.18) becomes 

^{vb z h)L\v(t)) + {yb y l h - y)L\v(t)) = {{zh)\y l h ) L ^Y(t)) + ((w - A)y ft , z h ) L 2 {vh{t)) + R h . 

Continue as in the proof of Theorem 5.11 to finally arrive at the analogue of (5.19) 
[z h A x {t)z h ]l + j \\y l h -y\\ 2 L 2 {m) dt<C(h 4 + sup || X 1 - /x|||^ (r(t)) ) J \\y h \\ 2 L2{rh(t)) dt . . . 

[0,T] * e[ °' T] [0,T] 

+ Ch 2 J ||y||^i (r(t)) + h 2 \\y h \\ 2 L2{rh{t)) dt, 

[0,T] 

Note that due to Lemma 5.5 

=0 since y%=Pfr((yo)i) 

\z h (0)A x (0)z h (0)\ = \(y l h (0)-y(0),z h ) L 2 (ro) \ < \(y h (0)-yi(0),z h } L2{rho) \+Ch 2 \\y \\ 2 L 2 {ro) . 

In view of Lemma 5.8 it remains to bound h 2 \\yh\\ 2 L2( j. h ^dt. Again thanks to Lemma 5.8 
we have 

f T 2 h 2 

\\yh\\ L 2(rh( t ))dt < C\\y \\ H i(r ) ■ 

But an inverse estimate, compare for example [CL91, Thm. 17.2], yields 1 1 2/o 1 1 J/ 1 (rg) — T^ll^o Hl 2 ( 
and because of the continuity of the lift (■); and of the L 2 -projection Pq the theorem follows. 
I 



6 Implicit Euler discretization 



In order to solve (3.8) we apply a vertical method of lines. The time discretization is carried 

2 

J L 2 (Y h y 



out by discontinuous Galerkin - implicit Euler discretization in L 2 2(Thy For N € N, consider 



an equidistant partition I n = (t n -i,t n ] of [0,T] with l<n<N,k=-^ and t n = kn. The 
trial space for the discontinuous Galerkin method (DGM) is the space of 'piecewise constant' 
functions 



W, 



t = [v G L 2 L2(rh) | VI < n < N : 3v n G V h (t n ) : v = 4 n:h v n on /„ } . 
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Note that in the following we will omit the operators <j>f h when dealing with functions w G W£. 

Also, to further simplify notation let a(t; ip, ip) = f V r hip ■ V T h<p + Xif)(pdr h (t) as well as 

i*(t) 

(■) -}n = (■, •) L 2 (v h (t n ))- W.l.o.g. we temporarily assume 

inf A( 7 ,t)>M + 2, (6.1) 
te[o,T\,ieTHt) 

with M = sup Te [ T ] || div r' i (T)^illL oo (r' i (r)) sucn that 

a{t; ip, ip) - M\\(p\\l 2(Th{t)) > M\ 2 H i {rh{t)) + IMli2 (r h (t)) 
for all t G [0,T], h > and all ip G Jf 1 ^*)). 

To motivate the DGM insert the Ansatz y k h (t) = Y!^=i ^hiv" 1 ^) G W k with V n G v h{tn) 
into (5.8). 

If one understands the time-derivative in (5.8) in a distributional sense, integration over time 
formally yields 

(y n - y n ' 1 ,f)n-i + J a(t; y n , <p) + (y n div Th V h , <p)L^(t)) dt = J (fh, ^) L 2 (T h (t)) d* , 

In In 

for smooth test functions ip. Instead, apply test functions ip G Wj} and use y n = (p n = to 
obtain 

J (y n div rh V h , <p) L 2 {rh{t)) dt = (y n , <p n ) n - {y n , >p n )n-i ■ 

In 

Finally, to arrive at a computable scheme, lump the Integral over a(t, -,-) and replace the 
right-hand side appropriately. For arbitrary parameters yft G Vh{0) and fh G L 2 L2( ^ h ^ we 
rewrite the scheme as 

y° f = y% , G W£ , 1 < n < N : 

(V],<p)n ~ (V^Mn-l + kOniy],^ = J (^Jh, ip) n dt , ^ 

In 

where y^, fh, and A are the same as in (5.8). For the approximation of the integral a n we 
assume a n (ijj, ip) = a(t n ; 4> t n h 4>, <Pt % h l P) + x n{4>, v), with a remainder 

\tn(ip,tp)\ < C t k\\4)\\ H i (Th(tn)) \\ip\\ H i {Th{tn)) . (6.3) 

One possible choice is t n = for 1 < n < N, but when it comes to approximating an adjoint 
equation such as (4.3) we will want to choose r more freely. In order to proof convergence of 
the scheme (6.2) in L^ 2 , rh . we make use of stability properties of the adjoint scheme 

z^ 1 = ZT , Vv? G W£ , 1 < n < N : 

{z n g , <p) n - {z n g + \v)n + ka n (<p, z n g ) = J (^ h g h , <p) n dt . ( 6 - 4 ) 

In 

with gh G L 2 L2 ^ Th y zt G Vh{T). In Section 7 it will be important that given snapshots 

{F h {t n )}n=i °f the surface (6.2) and (6.4) can be evaluated exactly for certain right-hand 
sides fh and gh, e.g. gh G W^. Let us introduce the mean value of a function y G L 2 L2 ^ Th ^ over 
an interval /„. 



26 



Lemma and Definition 6.1. Let <p s th denote the pullback operator associated to the flow 
as in Lemma 3.2 and let s G [0, T\. The mean value of a function y G ^^(pft) ^ s defined 
as y n (s) = \ J 4>'l h ydt for t G L n . Because 

in 

t, h ydt = J #, fc # ih j/dt = # ih J <t> r tth ydt, 

In In In 

y n does not depend on s G [0, T]. 

Similarly one could define the mean value of y G Wr if one were to investigate a horizontal 
method-of-lines approach. 

Now for yo = 0, zt = the schemes are adjoint in the sense 

N N 
k^2(fh,z g )n = k^2(g^,y f ) n , 

n=l n=l 

i.e. the discrete solution operators fh y/ and g^ h-> z 3 are adjoint as operators from 
{L 2 L 2^ vh y (-, into itself, where L^ 2 (ph) i s equipped with the scalar product 

N 

(f,g) h ,k = kJ2 ((^J)d^ h g))ndt. (6.5) 

n=l j 

Lemma 6.2. Let \\ ■ \\hk denote the norm induced by {-,-)hk- The norms \\ • \\ L 2 and 

L 2 (V h ) 



h t k on L L 2( V h} are equivalent and there holds 



(f,g)h,k- (f,g) L 2 9 , 
L 2 (r h ) 

Proof. The result follows from the identity 

N 



< Ck 



(f,g)L* . 

L 2 (r h ) 



J J fgdT h (t)dt = Y^J J (^J)(^ h g)Jl h dT h (t n )dt, 
ffli*(<) n=1 ^r l (i„) 



and J t tn being Lipschitz with j£ h = l. I 

Note also that for z G W£ , since i n = on L n , we can apply the mean value theorem to 
obtain for some t G L n 

\W zn \\h(rHt)) ~ ll*lnl = k \^ n d ™ r ^e n )Vh,z n ) L 2 {rh{@n)) \ < Mk\\z n \\l 2(rh{@n)) (6.6) 

with n G (t, i n )- Apply (6.6) to itself to obtain for some n G (0 ra , t n ) 

lllr n ll 2 II r ra ll 2 I <T A/f t /'ll ^™ll 2 J- /'ll y ra ll 2 1 1 7T. 1 1 2 \ \ 

\\\ Z \\L 2 (Th(t)) \\ z Ural ^ 1V1K [\\ z \\n-ry\\ z \\L 2 (r^(0 n )) W Z lira J J 

< Mk (\\z^l + Mk\\zX2 (rh( e n)) ) (6-7) 

< Mk(l + C L 2 (rh) Mk} \\z n \\ 2 n . 

A similar continuity assertion holds for the H l (T h (t))-norm, as shows the following lemma. 
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Lemma 6.3. Let y,z G , A G C(r /l (s) x [0,T]), and A = <^>* ft A. T/iere exisfe C > sttc/i 
i/tai /or ever?/ s G I n 

j *(s;<f> s t,hy,<t>t,h z ) dt - a(t;y,z)dt <Ck H s t , h y\\m(rh( s )) Wt,h z \\m(rH a )) dt > 
i.e. for z G W^ 1 we have 

ka(s;y n ,z n )- a(t;y,z)dt <Ck Ul h y\\m(rHs))\\z n \\m(rHs)) dt ■ 
j i n j i n 

In particular with A = 1 the estimates hold for a(t; (p, (p) = ll ( / 7 ll^i(r' i (t)) ' 

Proof. We abbreviate A(s,t) = D rh ( s ffi h (D rh ( s ffi h ) T J£ h . Since z n = we have 

I / a(s; #,hJ/, </>tV) d * - / a(t; y, z)dt| = . . . 

j 7 n J I n 

= 11 V r ^ t >(A(s, S )-A(M))v r ^ 
r h (s) 

The lemma follows from the fact that it linear on each T^(s) and globally Lipschitz in 
time, as by Lemma 5.3. I 

Let us formulate a crucial stability assertion for the adjoint scheme (6.4). 

Lemma 6.4. Let z G solve (6.4) with right-hand side g G L? L2 ^ h y For sufficiently small 
k > there exists C > 0, depending only on T, such that 

N N 

71=1 71=1 

Proof. Apply (6.4) to z n to obtain 

{z n _ z n+l^ z n )n + ^ = j (( ^ ^ dt . 

In 

This leads to 

2 (11*1* + \\* n+l ~ *ln ~ \\z n+1 fn) + ka n (z n , z n ) = I (<f% h g, z n ) n dt 

In 

< j ll^k Wt „))dt||^||n<^ (j U^OWndt) 
In V„ / 

Summing up and using (6.7) gives us 

E " *ln ~Mk(l + \c LHvh) Mk^ \\z n \\ 2 n + ka n (z n , z n ^j < ^ 
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n=l n=l ^ \ / y 



(6.8) 



< 



2M 1 



Now we test (6.4) with z n — z n+l to get 



\ v n _ _n+li|2 



+ - (a„(/\ z") + a n (z n+1 - - z n ) - a n (z n+1 , z n+1 )) =... 



(1/ 1 +-||z n -^ n+1 ||t 



= J(4>\: h 9^ n -z^) n dt<\(j WfyWl 
Summing up and using Lemma 6.3 on a as well as the estimate (6.3) on r we arrive at 
^.^O + ^E (\\z n+1 - z n \\l) 

n=m \ J 
1 k N 

< 2 k\\g\\ik+2 E a(tn-i;z n ,z n )-a(t n ;z n ,z n ) + x n ^(z n ,z n )-*n(z n ,z n ) 

n=m+l 
1 k N 

< 2 k \\9\\l, k + ^ E Ck (\\ zn Wm(r(t n )) + ll zn H|i(r(tn-i))) ■ 

n=m+l 

Combine with (6.8) to arrive at the lemma. I 

The following Lemma shows, that it is sufficient to estimate the approximation error at the 
points t n , 1 < n < N to prove convergence in L 2 L2 ^ Th y 

Lemma 6.5. Let r G H l ([0,T], V), V a separable Hilbert space, then there holds for r G I n 

\\r-r(T)\\ L 2 iIny) < k\\r'\\ L 2 {Iny) . 
In our situation this implies for r £ Hy h that 

1. k\\r(r) - f«||2 2(rh(r)) < Ck 2 f \\r\\ 2 L2{rh{t)) dt , 

In 

2. and f \\r(t) - f n \\ 2 L2[vh[t)) dt < Ck 2 f ||r||i 2(r , (t)) dt . 

in In 

Proof. For the fist assertion approximate r by n e V([0,T],V) such that rj — L >-' r as 
i — > oo. Use 

\\n-n(T)\\ L . (Iny) ={j jy^M v 6t\ < ( J k \\r' t (0)\\l d6 dt) < k\\ri\\ LHlny) , 
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and the fact that r G C([0,T],V), compare [LM68, Thm. 3.1]. Hence the first part of the 
lemma follows by passing to the limit. 

In our situation this implies, since (f)^ h r(t) G -ff 1 ([0, T], V^(r)) 



| f r i _ r(r M,2 



L3(rfc(T)) 



r(T) 



L 2 (r^(r)) 



dt 



In in 

This proves J., in order to get 2. integrate over I n . I 
We are now prepared to prove the main result of this section. 

Theorem 6.6. Let f G L 2 L2 ^ T y an d ^ Uh and y^^k solve (5.8) and (6.2), respectively, with 
yo ^ L 2 (Tq) and fh = fi- There exists a constant C > independent of h, k > and of f 
and y\ such that 



\\y h - y^hi Hrh) < Ck [\\y h \\Li 2(rh) + ll/ll^ 2(r) + WvoWl^) 

Proof. The proof is inspired by [SD05, Thm. 5.2], compare also [Vie07, Thm 1.2.5] and 
[MV08a, Thm 5.1]. Test (5.8) with <p\ n h ip, ip G Vu and integrate over I n to obtain 

(yh(t n ),<p)n- (yh(t n -i),<p)n-i+ J a(t; y h , ip) dt = J (fi,<p) L 2(rh(tydt. (6.9) 



[0,T] 



[0,T] 



Solve the adjoint equation (6.4) for z with both right-hand side and test function tp = g 



J M ~ vlkWl dt = (z n - z n+ \y n h ~ yl, k )n + ka n (y n h - yfo, 



(6.10) 



Subtract (6.9) from (6.2). Tested with z this yields 

(vtk - y h (t n ), z n ) n - (y^ 1 - y h (t n -i),z n ) n ^ + fca n « fc - yt z n ) = ... 

= j a(t;y hl z n )dt-ka n (yH,z n ) + k(fP,z n ) n - J {f h z n ) L 2 {rh{t)) dt 

In In 
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Let yh = J2n=i Vh^-in- Add (6.10) and sum up over 1 < n < N to get 



(fi,z)h,k- (fi,z) L 2 +J2 \\yh-yh,k\\ 2 n dt + / a{t;y h ,z n )dt-ka{t n -yl,z n ) = ... 

L(T) n=J / 

J-n J-n 

N 

= O + (Vh ~ yh(tn),Z n ) n ~ «^ - y h (tn-l), On-l " <^ + \ t h ~ <fc)„ 

n=l 

N 

= (yZk - Vh(tN), z n+1 )n - (y° h , k - y h (t ), ^>o + £ * n ) + (SJJ " z " " ^ +1 >™ 

n=l 

AT 

= £ fct^jft, z») + - y h (t n ), z" - z" +1 ) ft , 

n=l 

and finally, bringing to bear everything we have, i.e. the estimates from Lemma 6.3 for a, 
from Lemma 6.2 for the L 2 -norms, and the bound on r from (6.3), we arrive at 

£;/||Sfc-Vfc,*llndt< (kJ2M-y h (tn)\\l) (^I> n -* n+1 |lnl +"■ 
n=l j \ n=l / V n=l / 



J 2 II II r2 

L L2(r) V ( r) 
" «, ' 

<C|Mk* 



Hence using Lemma 6.4 on z we can divide by \\y~h — yh,k\\h,k- The Lemmas 6.2 and 6.3 allow 
us to estimate the involved norms, and because of the stability of the space discretization, 
compare Lemma 5.8, we can estimate the H 1 (T h (t))-teTm, to finally arrive at 



i 

N \ 2 



y h -y h M\Li 2(rh) <c^kJ2\\yh-yh(tn)\\ 2 n j +^11^+^^^) . (6.11) 



We now apply Lemma 6.5 [2.] to the error = yh,k — Vh and the averaged error = yh :k — Vh 
and sum up to obtain \\ek — ekWi 2 — Ck\\yh\\ L 2 . Combine with (6.11) and 6.5[i.] to 

L 2 (V h ) L 2 (V h ) 

estimate 

||e fe || L 2 <Ck\\y h \\ L 2 +||e fc || £ a <Ck(\\y h \\ L 2 + \\f\\ L 2 + \\ L 2 (r%) ) ■ 

L2(r h ) L2(V h ) L2{V h ) \ L 2 (rh) L 2 (r) y O'J 



With view of the stability assertions from (3.10) and Lemma 5.8 and together with Theorem 
5.11 we get the following Corollary. 

Corollary 6.7. In the situation of Theorem 6.6 let in addition A = fii and yo £ H 2 (To), 
and choose y\ as the piecewise linear interpolation of (yo)/- There exists a constant C > 
independent of h,k > and of f and yo such that 

\\y l h ,k - yhi 2(r) < c(h 2 + k) (bolero) + ll/ll^ 2(r) ) • 
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As addressed in Remark 5.12, it should be possible to relax the condition on yo into yo G 
i^ 1 (ro) using the L 2 (To)-projection or the L 2 (TQ)-projection Pq. 
But even in the case of low regularity we still get a uniform estimate. 

Corollary 6.8. In the situation of Theorem 6.6 let only yo G L 2 (Tq) hold while f = 0. Let 
further yft = /^((yo)/)- There exists a constant C > independent of h,k > and of yo such 
that 

h l h,k -v\\l\^ < c [ h + SU P l|A'-MlU<»(r(t)) + t J l|yo||L2(r )- 
L2 < r > \ te[o,T] « / 

Proof. Regarding Theorem 5.13 and 6.6 it remains to bound ||y/J|r2 . Like in the proof of 

L 2 (V h ) 

Theorem 5.13, using Lemma 5.8 and an inverse estimate, we arrive at the desired estimate. I 

In particular, for k > 0, choose k = nh? and A such that sup te j 0j T] 11^' ~~ A*llz,°° (r(t)) ^ C^ 1 to 
get an C(/i)-convergent scheme. 

Remark 6.9. Note that our freedom in the choice of r now allows us to finally drop the 
conditions on A and /x, respectively, in (5.7), (5.8), and (6.1). Let us assume we want to 
approximate the solution y of (5.7) with \i = 0, yo G ^(Tfi)), and / G L^^y Now 

yh,k G Wj! solves 

yl k = yt V<p€W£,l<n<N : 

(y T h,k,v)n - (<fc\¥>)n-i + k J y^(t n )Vlk ■ VrHt n )VdT h (t n ) = k{ft,<p) n , 

r>>(t n ) 

iff y h , m ,x = Zn=i e- Xtn vlkU n G W* A > solves 

Vh,k,\ = Vo, V<p€W£,l<n<N : 

{Vh,kx, <P)n ~ (<fe,A^)n-i + k J V rh(tn) ^ jfc • V rh(tn) (^ + Ay^dr^) + /ct n (y n , 

r h (t n ) 
= k(e- xt ^fj:^) n , 



with 



kt n ^,ip) = (e Xk -l-\k){^,cp) n + k(e Xk -1) y V r . (tn) V'-V r ^ tn) ^dr' l (t n ) 



r*(t„) 

Taking into account that \\e~ xt f(t) ~En=i e_A *" !/,»/(*) II L 2 2 } < fc ll/llL 2 2( } , we apply Corol- 
lary 6.7 to y h , m ,x and conclude ||y^ fc - y\\ L 2 < Ce XT (h 2 + k). 

L Z (T) 

7 Variational Discretization 



We now return to problem (P^) which has the advantage over (Ft), that its adjoint equation 
satisfies the regularity estimate (3.10). For (Ft) this is not the case iff yr G L 2 (T(T)) \ 
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J ff 1 (r(T)). In the spirit of [Hin05], let us approximate (P d ) by 
(P$) 



mm «eL2 2(rh) := U S d( u ) ~ (yd)l\\l,k + § 



s.t. a < u < b, 

with {T h (t)} tG[0tT] as in Section 5 and S% : (L\ 2 ^ h y (■, -) h ,k) -> (££2^, (■, -)h,fc), A ^ Vf is 
defined through the scheme 6.2 with A = and y$ = 0. We choose the scalar product (-, -}h,k 
defined in (6.5) in order to obtain a computable scheme to evaluate , namely (6.4) with 



N+l 



0. Given snapshots {T h (t n )}^ =1 , the product (-,-)h,k can be evaluated exactly for 
functions (fh G as well as for P\ a fi\(fh)- 

Let = |u G ^2( r ft) I o < ^ < & |- As in (4.2) the first order necessary optimality condi- 
tion for an optimum of (P^) is 

(au h + s¥ (S$u h - (y d )i),v-u h ) h , k >0, Vt> G J7 a d . (7.1) 

First note that as in the continuous case the (-, -^^-orthogonal projection onto U£ d coincides 
with the point- wise projection P[ a ^(v). Similar to 4.6 we get 

u h = P [aM (-^p5(«)) , pS(«) = sS* - • (7-2) 

Equation (7.2) is amenable to a semi-smooth Newton method that, while still being im- 
plementable, operates entirely in L? L2 ^ h y The implementation requires one to resolve the 

boundary between the inactive set T u (t n ) = {7 G T(t n ) | a < — ^p d (u)[y] < fr} and the ac- 
tive set A u {t n ) = T h (t n ) \ X u (t n ) for 1 < n < N. For details on the implementation see 
[HV10] and [HV11]. Note that in order to implement and according to (6.2) and (6.4) 
for right-hand sides in Wj}, again one only needs to know the snapshots {^ h (t n )}n=o- The 
solution of (Pj) converges towards that of (P d ) and the order of convergence is optimal in 
the sense that it is given by the order of convergence of and S% . 

Theorem 7.1 (Order of Convergence for (P^)). Let u G L 2 L2 ^ r y Uh G L 2 L2 ^ h ^ be the solutions 
of (P d ) and (P d )' respectively. Let C > 1. Then for sufficiently small h, k > there holds 



2a\\u l h - u\\ 2 L2 + \\y l h - y\\ 2 L2 < cfa ((O'sfOi " sj) (y - y d ),u- u{) 

{{■) l s } d\-)i - s d ) 



L h ( r) 



l Ihv) 



with y = S d u and y h = S^u h . 



Proof. The proof is a modification of the one from [HPUU09, Thm. 3.4], compare also [HV10]. 
Let Pyh (•) denote the (-, ■) /^-orthogonal projection onto U^ d . We have 

Ul = P [aM = P l*M (-^W') = P ^ d {~l Pd ^ 1 ) ■ 
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Since Uh € U^, from the characterization of P^h (•) it follows 



(--Pd(u)l ~ ui,Uh ~ ui) hjk < 0. 
On the other hand we can plug ui into (7.1) and get 

(au h + Pd(u h ),ui - u h ) hjk > . 
Adding these inequalities yields 

a \\ u l ~ u h\\l,k M (Pd( u h) ~ Pd{u)l) , ui ~ u h) h;k 

=(Pd( u h) ~ Sd*(v - Vd)i,ui - u h ) Kk + (S%*(y - y d )i -p d (u) h ui - u h ) h>k . 
The first addend is estimated via 

(Pd( u h) ~ (Sd)*(v ~ Vd)h ui - u h ) Kk = (y h - yi, S d ui - y h ) h)k 

= ~hh - yi\\l,k + (yh - vuSdUi - yi) h:k 
< -\hh - yi\\h,k + \\\ s d u i - yi\\l,k • 

This yields 

2a\\ui - u h f hk + \\y h - yi\\l,k ^ 2 (( S d*(-)i ~ (-)iS d )(y - z),ui - u h ) htk + \\S%ui - yi\\l,k- 

The claim follows for sufficiently small h,k > 0, using the equivalence of the involved norms 
stated in Lemma 6.2. I 



For the problem 

(Pr) < 



min «ei 2 r2/ „^ J ( n ) := U S t( u ) ~ (yT)i\? L 2 {T h {T)) + f NIl» 



s.t. a < u < b , 



L 2 (T h ) 



one can prove a similar result. Here the operator Sj< is the map fh — > yf{T), according to 
the scheme (6.2) with A = 0. 

Theorem 7.2 (Order of Convergence for (PjO). Let u G L 2 L2 ^ r y Uh € L 2 L2( ^ h ^ be the solutions 
of (Pt) and (Py), respectively. Let C > 1. Then for sufficiently small h,k > there holds 



IK - + IK " yfe ( r ( T)) < c ( 2 < ((-) 1 St(-)i -S* T )(y- y T ),u- u{) 



r 2 



+ 



((■) l sH-)i - St) 



L3(r(T)) 



to^/i y = Stu and yh = Sl^Uh . 
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Now as to the convergence of (j,-) l S^*(-)i — S^j, note that taking the adjoint does not com- 
mute with the discretization. Indeed, apply the scheme (6.2) to the adjoint equation (4.3), 
i.e. A = -( div r(tn) y)/ to get 

z^+i = o , Vv? G W£ , 1 < n < N : 



J (<j% h9h , Cf) n dt = <p) n ~ (Z% +I ,<p) n + ... 

J V r"(*„)^(i„) z 9 " - ( div r{t n )V)wz n dr fe (t„) + J {<p div rh(t) V h , z n ) L 2 {Th{t)) dt , 



In 

k 

r h (t n ) 



instead of (6.4). 

In the situation of (Pj) however, this discrepancy can be remedied by Lemma 5.6 which 
implies 

HO* - (■)'* \\c { l\ ( x 2 2i h ), IK-)' - (■Y\Wli 2 h ,L» ) < ch 2 , 

v L 2 (r)' L 2 (Y h )' v L 2 {T h ) L 2 (F)' 

and due to Lemma 6.2 which allows us to conclude 

11(0. - Onk^^x,.)^)), IK-)' - (■)«1lr(L ?JIfrfc „(^,^,(v>M)) * C ^ + ^ > ^ 



L 2 (r hyy^ L 2 (r hy 



if we interpret (•);, (•)' as operators into or on {L 2 L2 ^ h y (■, -}h,k), respectively. 
Hence we get the estimate 



< C(k + h 2 ) , 
in the C{L 2 L2 t T y ^^pp-operator norm. 

Now as to (Pt)> a h the results from section 5 and 6 remain valid under the time transform 
t' = T — t. As opposed to problem (P^), here it is easier to proof the convergence of S T * 
than that of S T itself. In order to discretize St* we choose A = — div r>»(t)^/i to approximate 
fil = -( div r(t) F)j and 



<p) = J((p div r fc (t ) Vh, ^ n ) L 2(rh(<)) dt _ div r "(t n ) Vh, z n ) L 2 



and apply Theorem 6.8 to end up with S T * (-)i — £tII£(l 2 (T(t)),l 2 2 ) — C(h+ |), where 
5f :z T 4zG^C ( L L2(rh)' (•> according to 



z T , V^ff,l<n<]V: 

(/,^ n -(/ +1 )¥) ) n + fc | v rh(tB) z n -v rh(tn) pdr h (i n ) = o. 

r^(t n ) 

Now in addition to (7.3) we have 

II ~~ (')' \\c(L 2 (r(T)),L 2 (r h (T))) ^ Ch 2 , 
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time=0.22 



time=0.44 



time=0.66 



time=1 




0.5 



-0.5 



Figure 1: Selected time snapshots of Uh computed for Example 8.1 on the Sphere after 4 
refinements. 



due to the inclusion (5.5). We conclude 

(•) 1 St(-)i — St 



C(Ll 2(r) ,LHT(T))) 



< C(h + 



the operator S T = Sf : {L^ hy (-, -) h , k ) -> L 2 (T h (T))), f h ^ y(T) being defined by the 
scheme 

y° = , \ftp € W% , 1 < n < JV : 

d/ n ,y)„- (l/"" 1 , ¥»)„-! + fc / V rh{tn) y n -V rh{tn) ipdT h (t n ) = k(fj;,v)n, 

T h (t n ) 

as shows summation over n. If is more regular, such as i/t £ -ff 1 (r(T)), then one might 
want to apply results from [DE11] that state /i 2 -convergence of the discretization S T , yet 
not in the C{L 2 L2{vy L?(T(T))) -norm. In order to to so, it remains to ensure the regularity 
assumptions of [DE11, Thm. 4.4] to be met by the optimal control u. 



8 Example 

Provided the results from [HIK03] and [Ulb03] hold on surfaces, Equation (7.2) is semi-smooth 
due to the smoothing properties of S%*, i-e. the stability ensured by Lemma 6.4. The lemma 
a priori holds only in the case A > 1, but can be extended for arbitrary A, /i by rescaling, see 
Remark 6.9. By Lemma 6.4 the operator cj) s h S% continuously maps {L 2 L2( ^. h y (•, -)h,k) m to 

L°°([0, T\,H l (T h (s))) C L p ([0, T], L' p (T h (s))) ~ L p ([0, T] x T h (s)) 
for every 2 < p < oo. This would imply semi-smoothness of the operator 

P[a,b] (-l<f>! lh ( P h d «,(•)))) : L 2 ([0,T] x T h (s)) L 2 ([0,T] x T h (s)) , 
compare [Ulb03], and thus of equation (7.2). 

We implemented a semi-smooth Newton Algorithm for (7.2), along the lines of [HV11]. 

Example 8.1 (High Regularity). Consider problem (P^) with a = 1, a = — i, b = |, 

T = 1, and T C M 3 the unit sphere. Let r(i) = ^qTq with i^(x,y, z) = (x , y , z / p 2 (t)) T and 
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6.78e-03 


2.15 


1.08e-01 


2.01 


1 


5.40e-02 




7.88e-01 




6 


3.15e-03 


2.09 


5.01e-02 


1.97 


2 


4.13e-02 


2.45 


5.32e-01 


0.86 


7 


1.72e-03 


2.03 


2.80e-02 


1.99 


3 


2.60e-02 


1.78 


3.78e-01 


1.79 


8 


7.92e-04 


2.02 


1.31e-02 


1.97 


4 


1.24e-02 


2.21 


1.82e-01 


1.97 













Table 1: L 2 -error, L°°-error and the corresponding EOCs for Example 8.1. 



sin(27Tt) 



pit) = e 2 ' . In coordinates (x, y, z) of M 3 let u = P^_i ij (z sin(27rt)) and yd = Vd + SdU 
with 



Vd = -a 



IT 

— sin 
2 



(2wt) - 2vr) 



.„ . sm(2irt)p(t) 

COs(27Tt) + > - 1 ^ 1 



X Z _|_ yZ _|_ 



x z + y 2 + p 2 z 2 



z . 



Then u solves (Prf). 



In order to compute the solution u/j of (P^) we construct triangulations of To from our macro- 
triangularion Rq, i.e. the cube whose nodes reside on To triangulated into 12 rectangular 
triangles. We generate Ri+i from Ri through longest edge refinement followed by projecting 
the inserted vertices onto IV 

Table 1 shows the relative error in the L 2 jj2^ h y\iovn\ and the relative -L°°-error 



ERR a 



^t,h( u h - u l)\\L°°([0,T]xr h (s)) 



HV ) i,/i' u 'llL oo ([0,T]xr' l (s)) 

as well as the corresponding experimental orders of convergence 

-l 



EOd = In 



ERRj 



ERR, 



In ( In 



i-q 



Hi 



Hi 



-q 



where H denotes the maximal edge length of Tq, see Table 2. Throughout this section we 
chose q = 2 for both EOC L 2 and EOCl°°, and the time step length is k = ^qH 2 . 
Figure 1 shows the solution of (P^) at different points in time. Note that the white line marks 
the border between active and inactive sets. On the active parts, the optimal control assumes 
the value a or b, respectively. 

Let us conclude with an example for (Pj,) with a desired state yr that just barely lies in 
L 2 (T(T)). In this situation we can only expect C(/i)-convergence. We consider the uncon- 
strained problem 

Example 8.2 (Low Regularity). Consider problem (Py) with a = 1, a = — oo, b = oo, T = 1 
and r(t) as in Example 8.1. Let yT ~ 



7TTT 



(x+y) 

Since we do not know the exact solution of Example 8.2, we estimate the relative error 



by ERR l L2 



U i+2\\L 2 



L 2 (rl +2)' 



N+2||l2 



i 2 (r»+2) 



where Ui denotes the solution of (P^) on 



the ith refinement {p i {t)}te[o,T} °f {^(^)}te[o,T]- The lift (■) is taken perpendicular to the 
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1 


2 


3 


4 


5 


6 


7 


ERR L 2 


0.1984 


0.0982 


0.0771 


0.0519 


0.0369 


0.0265 


0.0193 


0.0138 


E0C L 2 






1.6460 


1.5501 


1.3521 


1.0755 


0.9928 


0.9665 


H 


1.6330 


1.1547 


0.9194 


0.7654 


0.5333 


0.4099 


0.2769 


0.2085 



Table 2: L 2 -error and the corresponding EOC for Example 8.2. H is the maximal edge length 
of Tq (both examples). 

smooth surface T(t). Table 2 shows the estimated L 2 -errors and corresponding EOCs. We 
computed the L 2 (r fe (T))-projection P^Uti analytically. Otherwise the error introduced by 
the numerical integration of the non-smooth function yx would be dominant. It helps that 
all our triangulations resolve the plane {x + y = 0}. 
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